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ABSTRACT 

We present very high signal-to- noise ratio absorption-hne observations of CN and CH+ along 13 
lines of sight through diffuse molecular clouds. The data are examined to extract precise isotopologic 
ratios of ^^CN/^^CN and ^■^CH+/^^CH+ in order to assess predictions of diffuse cloud chemistry. Our 
results on ^■^CH+/^^CII+ confirm that this ratio does not deviate from the ambient ^^C/^^C ratio in 
local interstellar clouds, as expected if the formation of CH+ involves nonthermal processes. We find 
that ^^CN/^^CN, however, can be significantly fractionated away from the ambient value. The disper- 
sion in our sample of -'^^CN/^'^CN ratios is similar to that found in recent surveys of ^^CO/^'^CO. For 
sight lines where both ratios have been determined, the -'^^CN/^'^CN ratios are generally fractionated 
in the opposite sense compared to ^^CO/^^CO. Chemical fractionation in CO results from competition 
between selective photodissociation and isotopic charge exchange. An inverse relationship between 
i^CN/i^CN and i^CO/^^CO follows from the coexistence of CN and CO in diffuse cloud cores. How- 
ever, an isotopic charge exchange reaction with CN may mitigate the enhancements in ^^CN/^'^CN 
for lines of sight with low ^^CO/^'^CO ratios. For two sight lines with high values of ^^CO/^'^CO, 
our results indicate that about 50% of the carbon is locked up in CO, which is consistent with the 
notion that these sight lines probe molecular cloud envelopes where the transition from C^ to CO is 
expected to occur. An analysis of CN rotational excitation yields a weighted mean value for Tqi (^^CN) 
of 2.754 ±0.002 K, which implies an excess over the temperature of the cosmic microwave background 
of only 29 ± 3 mK. This modest excess eliminates the need for a local excitation mechanism beyond 
electron and neutral collisions. The rotational excitation temperatures in ^'^CN show no excess over 
the temperature of the CMB. 

Subject headings: ISM: abundances — ISM: clouds — ISM: molecules 
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1. INTRODUCTION 

The ^^C-to-^'^C isotopic ratio in the interstellar 
medium (ISM) is an important diagnostic for probing the 
history of nucleosynthesis and chemical enrichment in the 
Galaxy. The abundance of ^^C increases relative to ^^C 
over Galactic timescales because ^^C is a secondary prod- 
uct of stellar nucleosynthesis. Moreover, millimeter- wave 
emission studies of CO (Langer & Penzias 1990), II2CO 
(Henkel et al. 1982), and CN (Savage et al. 2002; Milam 
et al. 2005) in dense molecular clouds have demonstrated 
the existence of a Galactic gradient in the ^^C/^^C ra- 
tio, with increasingly lower values found closer to the 
Galactic center (see also the reviews by Wilson & Rood 
1994 and Wilson 1999). The increase in ^^C/^^C with 
Galactocentric distance presumably refiects the reduced 
rates of star formation, and subsequent stellar nucleosyn- 
thesis and mass loss, that are characteristic of the outer 
Galactic disk. The enrichment process, acting over the 
4.6 Gyr since the formation of the solar system, also offers 
a natural explanation for the reduction in the ^^C/^'^C 

^ Based in part on observations made with the Very Large 
Telescope of the European Southern Observatory, Paranal, 
Chile, under programs 065.1-0526, 071.C-0367, 071.C-0513, and 
076.C-0431. 
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ratio in the vicinity of the Sun. The terrestrial value of 
^^C/^'^C, which should approximate the interstellar value 
at the time of collapse of the solar nebula, is 89 (Lodders 
2003), while values close to 70 are typical of local inter- 
stellar clouds. Estimates for the present-day, local ISM 
value of ^^C/^'^C are usually obtained from absorption 
or emission studies of carbon-bearing molecules, which 
serve as proxies when direct measurements of the ambi- 
ent carbon isotopic ratio are lacking. Wilson (1999), for 
example, derived an average ^^C/^'^C ratio of 69 ±6 from 
millimeter- wave data on CO and II2CO emission for 13 
local sources. More recently, Milam et al. (2005) added 
millimeter observations of CN in dense clouds to those 
of CO and H2CO and found an average for the present- 
day, local ^^C/"'^^C ratio of 68 ± 15. In this paper, we 
shall follow Sheffer et al. (2007) and take the ambient 
12q^13q YdXio in the solar neighborhood to be 70 ±7 (see 
also Sheffer et al. 2002). 

When molecular proxies are used to evaluate the inter- 
stellar abundance ratio of ^^C to ^'^C, the values some- 
times show evidence of chemical fractionation. While 
this is generally not the case in dense molecular clouds, 
where the processes responsible for fractionation are less 
effective (see Milam et al. 2005), fractionation can have 
a significant impact in diffuse clouds. Detailed models 
predict that the fractionation processes in diffuse en- 
vironments will not affect each interstellar molecule in 
the same way. The molecular ion CII+ is not subject 
to chemical fractionation because it must be formed at 
high effective temperatures. The production of CII+ is 
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linked to the endothermic reaction C+ + H2 ^ CH+ + H 
(Elitzur & Watson 1978, 1980), which has an activa- 
tion energy of AE/kn = 4640 K. Since CH+ is asso- 
ciated with cold diffuse clouds, its formation requires 
that nonthermal processes, such as magnetohydrody- 
namic shocks or propagating Alfven waves, provide the 
additional heating. As a result, the ^^CH+/^''CH+ ra- 
tio is believed to be equilibrated and thus to be the best 
measure of the ambient carbon isotopic ratio in the dif- 
fuse ISM. Many observations (e.g., Centurion & Vladilo 
1991; Crane et al. 1991; Stahl & Wilson 1992; Centurion 
et al. 1995) have confirmed this expectation for CH+, 
revealing ^^CH+/^'^CH+ ratios very near 70 for local dif- 
fuse clouds. Centurion et al. (1995) give a weighted mean 
value of /^^CB.+ of 67±3 for five sight lines to stars 

within approximately 500 pc of the Sun. However, some 
authors (e.g., Hawkins & Jura 1987; Vladilo et al. 1993; 
Casassus et al. 2005; Stahl et al. 2008) have reported 
discrepant values of i2CH+/i-''CH+. Vladilo et al. (1993) 
found ratios of 126 ± 29 and 98 ± 19 toward HD 152235 
and HD 152424, respectively. Since both of these stars 
belong to the Sco OBI association (at d ~ 2 kpc), Vladilo 
et al. (1993) suggest that ^'^Cll+ /^^CB+ may vary on 
scales of ^1 kpc or smaller, depending on the distance to 
the clouds. Hawkins & Jura (1987) found ^^CR+/^^C}i+ 
ratios of 40 ± 9 and 41 ± 9 toward 20 Tau and 23 Tau, 
respectively. Because these stars are members of the 
nearby Pleiades cluster {d = 110 pc), such small val- 
ues of ^^CH+/^^CH+ seem to challenge the view that a 
ratio near 70 characterizes the solar neighborhood. 

Carbon-bearing molecules susceptible to fractionation, 
especially CO, may exhibit ^^C-to-^^C ratios that are 
either enhanced or reduced with respect to the ambi- 
ent value. Two competing processes are capable of al- 
tering the ^^CO/^^CO ratio in diffuse molecular gas. 
Selective photodissociation (SPD) favors -'^^CO since it 
is the more abundant isotopologue and is thus pro- 
tected to a greater extent via self shielding (e.g., van 
Dishoeck & Black 1988; Visser et al. 2009). The pro- 
cess is effective in the case of CO because the photodis- 
sociation of this molecule is governed by line absorp- 
tion. At lower gas kinetic temperatures, CO is influ- 
enced by the isotopic charge exchange (ICE) reaction 
13C+ + 12C0 ^ + 13C0 -I- A^;, where the differ- 

ence in zero-point energies AE/kB = 35 K (Watson et 
al. 1976) favors ^''CO. Because CO is the most abundant 
carbon-bearing molecule in the ISM, an enhancement in 
""^^CO, resulting from isotopic charge exchange, will de- 
plete the carbon reservoir of ^^C. Conversely, if ^^CO is 
selectively destroyed through photodissociation, then the 
carbon reservoir will be enhanced in ^^C. Any molecule 
arising from the remaining carbon in the reservoir should 
therefore possess a -'^^C-to--'^^C ratio that is fractionated 
in the opposite sense compared to the ratio in CO. Such 
behavior is expected for CN because this molecule co- 
exists with CO in diffuse molecular clouds (Pan et al. 
2005). Existing data on isotopologic ratios in CO, CN, 
and CH+ for the well-studied sight line to Oph offer 
some evidence that the chemical predictions are borne 
out. The diffuse gas in this direction has ^^CO/-^^CO — 
167 ± 15 (Lambert ct al. 1994), i^CN/^^CN = 47.3j:^-^ 
(Crane & Hegyi 1988), and ^^€11+ / ^^011+ = 67.5 ± 4.5 
(Crane et al. 1991). Thus, the ratios in CO and CN 
are fractionated in opposing directions, while the ratio 



in CH+ is consistent with the presumed ambient carbon 
isotopic ratio. 

In this investigation, we seek to measure ^^CN/^'^CN 
and ^^CH+/-'^'^CH+ ratios along lines of sight where 
^^CO/^^CO is either enhanced or reduced. Recent UV 
surveys of ^^CO/^^CO along diffuse and translucent sight 
lines (Sonnentruckcr ct al. 2007; Burgh ct al. 2007; Shcf- 
fer et al. 2007) have typically yielded ratios consistent 
with the average value of ^^C/ C for local interstellar 
clouds. Indeed, the weighted mean value of ^^CO/^'^CO 
for the 25 sight Unes studied by Sheffer et al. (2007) 
is 70 ± 2. However, there arc some notable exceptions. 
Enhanced ratios arc found, not only in the direction of 
C Oph. but also toward p Oph A and x Oph, where the 
respective ^^CO/^^CO ratios are 125 ± 23 and 117 ± 35 
(Federman et al. 2003), as well as toward ( Per, where 
i^CO/i^CO = 108±5 (Sheffer et al. 2007). Reduced ra- 
tios, with respect to the ambient value, are found in the 
directions of 20 Aql, where i^CO/^-'^CO = 50± 15 (Han- 
son et al. 1992), and HD 154368, where I^CO/^^CO = 
37 ± 8 (Sheffer ct al. 2007), among several others (see 
Sheffer et al. 2007). Sonnentrucker et al. (2007) find a 
low i^cO/i^CO ratio toward HD 73882 (25 ± 22), but 
the relative uncertainties are large. Our goal is to ob- 
tain i^CN/i-^CN and 12CH+/13CH+ ratios for as many 
of these sight lines as possible, so that, when our results 
are combined with the existing data on ^^CO/^'^CO, the 
suite of measurements can be used to test the predictions 
of chemical models for diffuse clouds. We reexamine the 
sight line to ( Oph as a check on our general method- 
ology. We also obtain ^^CH+/-^^CH+ ratios toward the 
Pleiades stars, 20 Tau and 23 Tau, to investigate poten- 
tial scatter in the carbon isotopic ratio within the solar 
neighborhood. 

Since high-quality data are needed to detect the weak 
features associated with '^'^CN and '^'^CH+, the observa- 
tions examined here also allow precise determinations of 
rotational excitation temperatures in CN. In some situa- 
tions, these measurements can then be used to constrain 
the electron density, an important parameter for mod- 
eling the physical conditions within a given cloud (e.g.. 
Black & van Dishoeck 1991). It is widely understood 
that the CN molecule is maintained in radiative equi- 
librium with the cosmic microwave background (CMB) 
in interstellar space (see the review by Thaddeus 1972). 
The CMB is the primary source of radiation in the uni- 
verse at 2.64 mm and 1.32 mm, the wavelengths of the 
two lowest rotational transitions in CN. As a result, CN 
excitation temperatures reflect the tcmperatiuc of the 
CMB at these wavelengths, in the; abscincc; of any local 
sources of excitation. Observed excitation temperatures 
in CN, in fact, do exhibit an excess over the temperature 
of the CMB, as derived from the Far Infrared Absolute 
Spectrophotometer (FIRAS) onboard the COBE satel- 
lite (TcMB = 2.725 ± 0.002 K; Mather et al. 1999; see 
also Fixsen 2009), but this excess is small (< 0.1 K; e.g., 
Palazzi et al. 1992). Since electron impact should domi- 
nate any local contribution to CN excitation (Thaddeus 
1972), an observed excess provides an estimate for the 
density of electrons in the portion of the cloud traced by 
CN. 

In order to investigate isotopologic ratios in CN and 
CH+ and CN rotational excitation in diffuse molecular 
clouds, we examine high-resolution, very high signal-to- 
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noise ratio observations of optical absorption lines arising 
from electronic transitions within the CN B ^S"*"— X ^E"*" 
and CH"*" A ^11— X ^S"*" systems. The observations and 
data reduction procedures are described in § 2. In § 3, 
we provide detailed information concerning the profile 
synthesis routine, with which we derive our final column 
densities and isotopologic ratios. The analysis and dis- 
cussion of our results on ^^CIl+ /^^C}i+ and i^CN/^^CN 
ratios appears in §§ 4.1 and 4.2, respectively. In § 4.3, 
we explore the relationship between ^^CN/^'^CN and 
^^CO/^^CO in an effort to evaluate the effects of chem- 
ical fractionation in diffuse molecular gas. The topic of 
CN rotational excitation is examined in § 5, and our main 
findings are summarized in § 6. An appendix gives re- 
sults on weak Ca I and Ca II absorption toward a Leo 
and a Vir, which were observed as unreddened standard 
stars to aid in the reduction of the McDonald Observa- 
tory data (see § 2.1). 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. McDonald Observatory Spectra 

Principal data for this project were acquired with the 
Harlan J. Smith 2.7 m telescope at McDonald Observa- 
tory using the TuU (2dcoude) spectrograph (Tull et al. 
1995) in its high-resolution mode (cs21). The observa- 
tions were carried out over the course of four observing 
runs between 2007 January and 2008 May. Two winter 
runs (2007 January and 2007 December/2008 January) 
provided data on 20 Tan, 23 Tau, and ( Per, while spec- 
tra of p Oph A, C Oph, and 20 Aql were acquired dur- 
ing two summer observing sessions (2007 June and 2008 
May). For all of the observations, the crossed-dispersed 
echelle spectrometer was configured with the 79 gr mm^^ 
grating (El), the 145 ^m slit (Slit 2), and a 2048 x 2048 
CCD (TK3). Order 56 of the El grating spectrum was 
centered at 4077 A, enabling both the CN B-X (0, 0) 
lines near 3874 A and the CH+ A-X (0, 0) transition 
at 4232 A to be recorded with a single exposure0. The 
grating tilt was adjusted slightly each night so that the 
interstellar features would not always fall on the same 
portion of the CCD. This helps to mitigate the effect of 
fixed pattern noise and reduces the chance that an in- 
strumental glitch at a specific location on the CCD will 
negatively impact the observations. 

The detection of absorption from the ^^C-bearing iso- 
topologues of CN and CII+ requires spectra with very 
high signal-to- noise ratios (S/N > 1000 per resolution el- 
ement). To achieve this, the program stars were observed 
for several hours each night. Individual stellar exposures, 
taken in succession, were limited to 30 minutes to mini- 
mize the number of cosmic ray hits in each integration. 
A bright, unreddened star was observed nightly to the 
same signal level as the program stars to confirm that 
no instrumental defects were present at the locations of 
interstellar features. For both of the winter observing 
runs, a Leo (Regulus; B7 V; B = 1.24) served as the 
unreddened standard star. For the summer observations, 
either a Vir (Spica; Bl III-IV-I-; B = 0.91) or, on one 
occasion, a Lyr (Vega; AO Y; B — 0.03) was observed for 
this purpose. Calibration exposures for dark current and 

^ This setting also provides data on the CH A—X (0, 0) transi- 
tion at 4300 A and on the Ca I A4226 and Ca II A3933 hnes. 



bias correction, fiat fielding, and wavelength assignment 
were obtained during each of the observing runs. Four 
30-minute dark frames were acquired on the first night 
of a run, while a series of ten biases and ten fiats were 
taken each night. For accurate wavelength calibration, 
Th-Ar comparison spectra were recorded throughout the 
night at intervals of 2— 3 hr. 

Reduction of the raw McDonald data employed stan- 
dard routines within the Image Reduction and Analysis 
Facility (IRAF) environment. First, the overscan region 
in each of the raw images was fitted with a low-order 
polynomial and the excess counts were removed from the 
exposure. An average bias was created for each night and 
was used to correct the darks, flats, stellar exposures, and 
comparison lamp frames. No dark correction was applied 
to any of the observations because the level of dark cur- 
rent was always found to be insignificant after the bias 
was removed. Cosmic rays were eliminated from the stel- 
lar and comparison lamp exposures by setting the respec- 
tive thresholds for detection to 50 and 25 times the mean 
of the surrounding pixels. Any cosmic rays present in in- 
dividual flat lamp exposures were effectively removed by 
taking the median of all flats for a given night. Scattered 
light was modeled in the dispersion and cross-dispersion 
directions and subtracted from the stellar exposures and 
from the median flat. The flat was then normalized to 
unity and divided into the stellar and comparison lamp 
frames to account for (flrst-ordcr) pixel-to-pixel varia- 
tions across the CCD. One-dimensional spectra were ex- 
tracted from the processed images by summing across 
the width of each order without weighting the individual 
pixels. 

At this point, the extracted spectra of the program 
stars were examined in pixel space to determine the 
pixel coordinates of the interstellar absorption lines. The 
combined spectrum of the unreddened star from each 
night was then inspected at these positions to deter- 
mine whether any instrumental defects remained after 
flat flelding. No obvious glitches were found in the vicin- 
ity of CN or CH+. However, low level fixed pattern noise 
was identified in some of the stellar spectra, particularly 
those with the highest S/N. To correct for these second- 
order fiuctuations, the affected stellar exposures were di- 
vided by the normalized spectrum of the unreddened star 
for that night. The same technique, when applied to ex- 
posures with lower S /N, did not significantly improve the 
quality of the data and so was not adopted for these expo- 
sures. Next, stellar spectra were calibrated in wavelength 
after identifying emission lines in the Th-Ar comparison 
spectra, typically 5 per order. The dispersion solution 
applied to a given stellar exposure resulted from an in- 
terpolation between the solutions derived from the Th- 
Ar lamp exposures that preceded and followed the stellar 
exposure. Finally, the calibrated spectra were shifted to 
the reference frame of the local standard of rest (LSR) . 

All exposures of a given target from the multiple ob- 
serving runs were coadded to maximize the S/N in the 
final spectrum. Before inclusion in the sum, each ex- 
posure was carefully examined, speciflcally near the ex- 
pected positions of the ^'^CN and ^'^CH+ absorption fea- 
tures, to determine if any weak cosmic rays had survived 
the reduction process. When cosmic rays were identi- 
fled, they were removed manually if they occurred far 
enough away from the interstellar features. If the cosmic 
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Figure 1. Possible continuum fits {dotted lines) to the region 
surrounding the CH+ lines toward 20 Tau. The strong interstellar 
^■^CH+ line is superimposed on the blue wing of a photospheric 
Fe II line, complicating the placement of the continuum. The po- 
sition of the ^"^011+ line is indicated by a tick mark. 



ray affected an interstellar feature directly, then the ex- 
posure was removed from the sum. Because the CN and 
CH+ lines reside in different echelle orders and the orders 
were coadded separately, an exposure removed from the 
sum in one order would not necessarily be removed in the 
other. Thus, the total exposure times for the two spectral 
regions may differ somewhat (see Table 1). Further com- 
plicating the coaddition process was a slight change in the 
resolving power of the spectrograph from one run to the 
next. The measured widths of thorium emission lines in 
the Th-Ar comparison spectra indicate a resolving power 
of i? = 173,000 (Av = 1.7 km s'^) was achieved dur- 
ing our run in 2007 January. This decreased to 150,000 
(2.0 km s~^) for the run in June of that year, but re- 
turned to 160,000 (1.9 km s'^) and 168,000 (1.8 km s'^) 
for the runs in 2007 December and 2008 May, respec- 
tively. When combining exposures from multiple observ- 
ing runs, the data were sampled at the dispersion of the 
lowest resolution spectra contributing to the sum. 

Final coadded spectra were normalized to the contin- 
uum by fitting low-order polynomials to regions free of 
interstellar absorption within a narrow spectral window 
(no more than 2 A wide) surrounding a given interstellar 
line. The immediate vicinity of the CN and CH+ lines 
in our McDonald spectra arc free of contamination by 
either telluric or instrumental features, making contin- 
uum normalization straightforward. In only one case did 
the placement of the continuum present any significant 
difficulties. The interstellar ^^CH+ line toward 20 Tau is 
superimposed onto a narrow photospheric Fe II feature 
at 4233 A. This situation is unique among the McDon- 
ald targets to 20 Tau, which has a relatively late spectral 
type and a rather slow rotation rate {v sin i ^ 39 km s~^; 
Hoffleit & Jaschek 1982). Figure 1 shows a range of pos- 
sible continuum fits to the regions surrounding the CH+ 
lines toward 20 Tau. The small variation in the contin- 
uum above the ^'^CH+ feature at 4232.4 A has the great- 
est impact on the resulting ^^CH+/^^CH+ ratio, but this 
variation is less than 10% of the ^•^CH+ line strength, 
which is well within the 20% observational errors associ- 
ated with the detection of this weak feature. 



2.2. VLT/UVES Archival Spectra 

To supplement our high signal-to-noise McDonald ob- 
servations of bright targets, optical spectra of seven 
fainter stars were obtained from the European Southern 
Observatory (ESO) Science Archive Facility. These spec- 
tra were acquired by the Ultraviolet and Visual Echelle 
Spectrograph (UVES) of the Very Large Telescope 
(VLT) at Cerro Paranal, Chile primarily during two ob- 
serving programs designed to measure -"^^011+ /-^^CH+ 
ratios (see Casassus ct al. 2005; Stahl ct al. 2008). 
The primary UVES data (from programs 071.C-0367 and 
076.C-0431; S. Casassus, PI) have S/N ratios that are 
comparable to those of our McDonald spectra. These 
observations employed the central wavelength setting at 
4370 A, which allows continuous spectral coverage from 
3750-5000 A, a range that includes the CN (0, 0), 
CH+ (0, 0), and CH+ (1, 0) lines. Other VLT/UVES 
datasets were obtained (from program 065.1-0526; E. 
Roueff, PI) that were acquired using the wavelength 
setting at 3460 A. This setting covers the region from 
3050-3870 A and permits detection of the CN (1, 0) 
lines near 3579 A. The observations at 3460 A typically 
have lower S/N than those at 4370 A. For one sight line 
(HD 210121), observations at 3460 A and 4370 A were 
provided by program 065.1-0526 and additional data ac- 
quired with the wavelength setting at 3900 A were ob- 
tained from program 071.C-0513 (M. Andre, PI). 

The weaker CN and CH+ (1, 0) lines provide a check 
on the column densities derived from the stronger (0, 0) 
lines, which is especially important when the (0, 0) lines 
become optically thick. The seven UVES sight lines cho- 
sen for this study (see Table 1) are known from previous 
observations (e.g., Gredel et al. 1991; Palazzi et al. 1992) 
to have strong lines of interstellar CN and determining 
-'^^CN/^'^CN ratios will be our focus. We do not rederive 
12CH+/13CH+ ratios for the UVES sight lines as there 
are significant uncertainties associated with such deter- 
minations from these data (see § 3). While the principal 
reason for incorporating the UVES spectra into our anal- 
ysis was to expand our sample to include more heavily- 
reddened sight lines, the VLT observing programs also 
provided us with additional data on some of our McDon- 
ald targets. Both 23 Tau and C Oph were observed at 
4370 A under programs 076.C-0431 and 071.C-0367, re- 
spectively, and observations of p Oph A and C Oph at 
3460 A arc available from program 065.1-0526. These 
supplementary data are useful for making comparisons 
with our McDonald results. 

All VLT/UVES datasets were reduced with the UVES 
pipeline software in optimal extraction mode unless the 
observations employed an image sheer in which case the 
average extraction method had to be adopted. Mas- 
ter bias and master flat frames were created for each 
night for which we obtained science observations. Th- 
Ar lamp spectra acquired contemporaneously with the 
science data were used for wavelength calibration. Af- 
ter extraction, all echelle orders from a given exposure 
were merged and the resulting spectra were shifted to the 
LSR frame of reference. Finally, all individual exposures 
of the same target obtained with the same instrumen- 
tal setup were coadded. The spectral resolution of the 
UVES data varies somewhat for the different observing 
programs due to changes in slit width and CCD binning. 
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Figure 2. Coadded UVES spectrum of CN B-X (0, 0) toward 
HD 169454 showing the adopted continuum fit {dotted line). Each 
line is labeled and the positions of the R(0), R{^), and P(l) lines 
of ^'^CN are indicated by tick marks. This is the first reported 
detection of the ^^CN P(l) line in the interstellar medium. Note 
also that absorption from ^'^CN R{0) is partially blended with the 
much stronger -"^^CN R(0) line. 

The primary data at 4370 A have a velocity resolution of 
Av = 3.5 km s~^, corresponding to a resolving power of 
R = 85, 000. The observations at 3460 A employed either 
1x1 or 1x2 binning and have resolutions of 3.8 km s~^ 
{R = 79,000) or 4.2 km s^^ (i? = 71,000), respectively. 
For HD 210121, the observations at 4370 A were com- 
bined with those at 3900 A in the overlapping region, 
and the resulting spectrum has a resolution correspond- 
ing to i? = 80,000. 

The coadded UVES spectra in the vicinity of CN were 
normalized to the continuum in much the same way as 
were the McDonald data except that a larger spectral 
window was used (typically 5 A) that included all of the 
hues of the CN (0, 0) or (1, 0) band. The R{0) and 
i?(l) lines in these data are separated by very little con- 
tinuum due to the lower resolution. The larger spectral 
windows, in addition to a non- negligible amount of resid- 
ual fringing in the extracted spectra, meant that higher- 
order polynomials had to be applied during the normal- 
ization process. Moreover, the ^^CN and ^■^CN i?(0) lines 
of the (0, 0) band, which are separated by approximately 
13.5 km s~^, are partially blended in the UVES spectra, 
again a result of the coarse resolution but also due to the 
strength of CN absorption along these lines of sight. Ab- 
sorption from "'^'^CN is thus found on the redward wing 
of the ^^CN line (as shown in Figure 2), adding some 
uncertainty to the placement of the continuum in this 
region. The situation is not as problematic as in the case 
of CII+ from these data, however, because the CN lines 
are intrinsically narrow and exhibit mostly single velocity 
components. 

2.3. Stellar Sample and Observational Results 

The total exposure times and signal-to-noise ratios per 
resolution element achieved for each of the McDonald 
and UVES targets are shown in Table 1 along with the 
relevant stellar data, consisting of the spectral type, B 
magnitude, E{B — V) color excess. Galactic / and b co- 
ordinates, and distance. We list the Hipparcos parallax 
distance (Ferryman et al. 1997) if the measurement is 
significant at the 4-a level or greater. Otherwise, the dis- 



tance we derive from the method of spectroscopic paral- 
lax is given. The spectral types in Table 1 are from the 
same references that provided values of E{B — V), al- 
though, for HD 73882, HD 154368, and HD 170740, the 
spectral types were actually listed in a previous paper 
by the same group (Rachford et al. 2002). From these 
references, we also obtained the V magnitude and total 
visual extinction (Ay) for each star without a precise 
Hipparcos measurement so that spectroscopic distances 
could be derived using a self-consistent set of photomet- 
ric data. As the distances in Table 1 indicate, our selec- 
tion of interstellar sight lines samples almost exclusively 
the local ISM. Only four stars in the sample (HD 73882, 
HD 152236, HD 154368, and HD 169454) lie farther than 
500 pc from the Sun. Accordingly, these stars also ex- 
hibit the highest degrees of interstellar reddening, with 
the largest value of E{B — V) found in the direction of 
HD 169454. Even along the more extended sight lines, 
however, the clouds responsible for CN and CH+ absorp- 
tion are likely to be nearby. This can be deduced from 
the velocities of the molecular components, which corre- 
spond to the dominant components in species like Ca II. 
The more distant stars have rather more complicated 
profiles of Ca II A3933, but the dominant components 
are always found near i^lsr = km s~^, implying that 
the absorption originates in local gas. 

The S/N ratios in Table 1 result in uncertainties in 
equivalent width (Wx) of 0.06 mA, on average, for the CN 
and CH+ lines in our McDonald spectra. For the primary 
UVES data, the average uncertainty in Wa(CN) is 0.05 
mA. Despite the small uncertainties, we do not detect 
absorption from ^'^CH+ toward ( Per or 20 Aql, which 
is not unexpected considering the weak ^^CH+ lines in 
these directions. The 3-tT upper limits on the equivalent 
widths of 13CH+ are < 0.18 mA for C Per and < 0.24 
mA for 20 Aql. We also do not detect any CN absorp- 
tion from either isotopologue in the directions of 20 Tau 
and 23 Tau, but obtain upper limits on Wa(^^CN) of 
< 0.25 mA and < 0.13 mA for the two sight lines, re- 
spectively. The equivalent widths for detected lines are 
given in Table 2 for CH+ A4232 and Table 3 for CN 
A3874. We list the fitted equivalent widths that result 
from our profile synthesis analysis (see § 3). Uncertain- 
ties in W\ were calculated as the product of the rms 
variations in the continuum and the full width at half 
maximum (FWHM) in the linc0. Note the quite strong 
detections of ^^CN toward the UVES targets. While not 
hsted in Table 3, the weak R{2) and P(2) lines of ^^CN 
are detected toward most of the stars in our sample with 
CN absorption. In our McDonald spectra, however, these 
features are only marginal. For C Oph and HD 152236, 

^ While we do not formally include continuum placement errors 
when calculating uncertainties in Wx, the uncertainties derived 
using our method are more conservative than those based on sta- 
tistical errors (adopting the formula given in Jenkins et al. 1973) 
and our estimates for continuum placement errors added in quadra- 
ture. As an example, the statistical error in the equivalent width 
of the line toward 20 Tau is awrms = (AA)(Afp'^)(rms) = 

0.0496 mA, where AA = 0.00802 A is the wavelength disper- 
sion, A'^p = 21 is the number of pixels sampling the line pro- 
file, and rms = 0.00135 is the root mean square of the noise in 
the continuum. The standard deviation in the measured value 
of Wa(^'^CH+) for a range of possible continuum fits (see Fig- 
ure 1) is (7 vi/jCont — 0.0229 mA. Thus, the total error would be 
o-w = (o-wrms + o'w'cont)^''^ = 0.0546 mA, while our method 
gives aw = '(FWHM)(rms) = 0.0730 mA. 



6 



RITCHEY, FEDERMAN, & LAMBERT 



20 Tau 
1.2 - 




I , , , I , ^ I 

4232.4 4232.6 4232.8 

Wavelength (A) 



Figure 3. Simultaneous profile synthesis fit to the CH+ and 
^•^CH+ R(0) lines toward 20 Tau scaled to show the stronger {up- 
per panel) and weaker (lower panel) isotopologuc. The synthetic 
profile is shown as a solid line passing through data points that 
represent the observed spectrum. Residuals are plotted above the 
fit. Individual velocity components are indicated by tick marks. 

the lower signal-to-noise ratios near 3874 A, compared to 
other stars observed with the same instrument, prevent 
any detection of i?(2) or P(2). Remarkably, four of the 
UVES targets exhibit fairly significant absorption from 
i^CN and in the spectrum of HD 169454, the ^^CN 
P(l) line is detected with a significance of 8 cr (see Fig- 
ure 2). Also shown in the tables are previous results from 
the literature on CH+ and CN toward our targets, where 
we have included only those studies that give equivalent 
widths for the weaker isotopologue or for the R{1) and 
P(l) lines in the case of CN. Studies that examined only 
the strongest line in CH+ or CN are omitted. 

3. PROFILE SYNTHESIS 

The CN and CH+ absorption profiles were synthe- 
sized with the rms-minimizing code ISMOD (Y. ShefFcr, 
unpublished) that treats the velocities (v), Doppler b- 
values, and column densities (A^) of the absorption com- 
ponents as free parameters. ISMOD assumes a Voigt 
profile function for each absorption component and con- 
volves the intrinsic profile with an instrumental profile, 
represented by a Gaussian with a width determined by 
the resolving power of the instrument (see Black & van 
Dishoeck 1988). The i?(0) lines of the (0, 0) bands in CN 
and CH+ were used to derive the N{^'^CN) /N(}^CN) and 
Ar(i2CH+)/iV(i3CH+) ratios. Both isotopologues of CN 
and of CH"'" were fitted simultaneously by requiring the 
component structure of the two profiles to be identical. 
Thus, the 12-to-13 column density ratio and the isotope 
shift between the two lines were additional free parame- 
ters. 

Figures 3 through 10 present the CH+ and CN R{0) fits 
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Figure 4. Same as Figure 3 except for the ^■^CH+ and ^"^013+ 
R{0) lines toward 23 Tau. 
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Figure 5. Same as Figure 3 except for the ^•^CH+ and ^^"^013+ 
R{0) lines toward p Oph A. 

for the data acquired at McDonald Observatory. When 

fitting the CH+ lines toward 20 Tau and 23 Tau, we used 
our previous results on the component structure in these 
directions (Ritchey et al. 2006) as the initial input to 
the profile synthesis routine. Likewise, the initial values 
for the component structure of CH+ toward p Oph A 



CN AND CH+ IN DIFFUSE MOLECULAR CLOUDS 



7 



f Oph 
1.2 - 




4232.2 4232.4 4232.6 4232.8 

Wavelength (A) 

Figure 6. Same as Figure 3 except for tlie ^^CH'^ and ^"^011+ 
R{0) lines toward C Oph. 



were obtained from Pan et al. (2004). The input values 

for the CN i?(G) hnc toward ( Oph were taken from the 
ultra-high resolution (i? = 600, 000) study by Lambert et 
al. (1990), who modeled CN absorption along this sight 
line with two Gaussian components separated by 1.2 km 
s~^. Because our study has coarser resolution and the 
component parameters in this case were already precisely 
determined, the velocity difference between the two com- 
ponents, their relative strengths, and their 6- values were 
held fixed in our analysis. Lambert et al. (1990) also 
examined CH+ toward ( Oph, but modeled the line with 
a single Gaussian component. Absorption profiles with 
single components do not require precise input parame- 
ters in ISMOD. Thus, for the CH+ hnc toward ( Oph, as 
well as for the CN i?(0) lines toward the other McDon- 
ald targets, input parameters were determined by direct 
examination of the spectra within IRAF. The velocities 
and column densities of the interstellar absorption com- 
ponents in our profile synthesis fits were computed using 
the wavelengths and oscillator strengths (/-values) pre- 
sented in Table 4. 

Among the optical tracers of diffuse molecular gas, the 
CN molecule typically probes the densest portion of a dif- 
fuse cloud (Pan et al. 2005). As a consequence, CN line 
widths arc intrinsically small, corresponding to 6-values 
that are less than 1.0 km s^^ in most cases. With such 
narrow lines, the strongest transitions in CN can become 
optically thick even at moderate column densities. While 
this is not a major concern for our McDonald sight lines, 
all of which have optical depths (r) at line center that 
are less than 1.0, it is an important consideration for the 
UVES observations, where the average value of t for the 
(0, 0) R{0) line is 4.45. When large optical depths are 
present, a small change in 6(CN) will have a dramatic 
impact on the ^^CN column density obtained from the 




3874.5 3874.8 3874.7 3874.8 3874.9 3875.0 
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Figure 7. Simultaneous profile synthesis fit to the ^'^CN and 
^^CN R{0) lines toward f Per scaled to show the stronger {upper 
panel) and weaker (lower panel) isotopologue. Plotting symbols 
are the same as in Figure 3. 
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Figure 8. Same as Figure 7 except for the ^^CN and ^^CN R{0) 
lines toward p Oph A. 

strong i?(0) line, and the resulting ^^CN/^'^CN ratio wiU 
be similarly affected. Moreover, because the UVES spec- 
tra are sampled at moderate resolution (~3.5 km s~^), 
the CN lines are dominated by instrumental rather than 
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Figure 9. Same as Figure 7 except for tlie ^^CN and ^^CN R{0) 
lines toward C Oph. 
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Figure 10. Same as Figure 7 except for tiie ^^CN and ^^CN R{0) 
lines toward 20 Aql. 



intrinsic broadening, making direct evaluation of the 6- 
valucs more difficult. To minimize these uncertainties, 
6-values for optically thick lines were derived using the 
well-known doublet ratio method, which has often been 
employed in this context (e.g., Meyer & Jura 1985; Gre- 
del et al. 1991; Roth & Meyer 1995). The method is 
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Figure 11. Column density as a function of b- value for the four 
CN transitions used to derive 6 in the main component toward 
HD 154368 via the doublet ratio method (sec text). The intersec- 
tion of each pair of curves, which gives the 6- value corresponding to 
the measured doublet ratio, is marked by a solid point. The I-ct er- 
rors on b thus derived arc shown as dashed lines for the (0, 0)/(l, 0) 
R{0) pair and dot-dashed lines for the (0, 0) R(1)/P{1) pair. 



predicated on the fact that two lines arising from the 
same state must yield identical column densities for some 
value of b. In the; case; of CN. one can use either the (0, 0) 
i?(l) and P(l) lines or the (0, 0) and (1, 0) E(0) lines. 
Having observations of both pairs is advantageous as it 
allows independent checks on consistency. 

In the optically thin limit, the equivalent width of an 
absorption line is proportional to the column density 
such that Wx oc X^fN (see Spitzer 1978). Thus, for a 
pair of optically thin lines that must yield identical values 
of N, the equivalent width ratio is given by VKa,i/Wa,2 = 
(Ai/A2)^/i//2- Specifically, from the molecular line pa- 
rameters in Table 4, the (0, 0) i?(l)/P(l) equivalent 
width ratio for optically thin absorption is equal to 2.00, 
while the (0, 0)/(l, 0) i?(0) ratio is equal to 13.35. The 
measured ratio will be lower if the optical depth of the 
stronger line is significant. To derive 6-values using this 
approach, the CN equivalent widths were determined by 
fitting simple Gaussians to the observed profiles within 
IRAF. Column densities based on the measured equiva- 
lent widths for a given pair of lines were then interpolated 
from a set of curves of growth calculated on the interval 
0.30 <b< 0.99 km s'^ with a step size of 0.01 km s'^ 
The result of this procedure was a pair of curves that give 
TV as a function of b for the two transitions, corresponding 
to the measured values of Wx. By plotting these curves 
and finding the intersection in the; N-b plane, one finds 
the unique 6- value that yields identical column densities 
for the pair of lines (see Figure 11). 

Table 5 gives the observed CN (0, 0) i?(l) /P(l) equiva- 
lent width ratios, and (0, 0)/(l, 0) R{0) ratios if possible, 
for all of the stars in our sample where CN is detected. 
For the main CN component toward HD 152236, both 
ratios are consistent with the optically thin limit. The 
iJ(l)/P(l) ratio of 1.91 ± 0.10 nominally gives a 6- value 
of 0.39 km s~^, but with an upper bound at infinity. An 
infinite 6-valuc in this context means that the lines fall 
on the linear portion of the curve of growth and cannot 
be used to constrain the Doppler parameter. This is ev- 
idently the case for the CN lines toward our McDonald 
targets, which yield no definitive information on b. It is 



CN AND CH+ IN DIFFUSE MOLECULAR CLOUDS 



9 



worth noting, however, that the (0, 0)/(l, 0) -R(O) equiva- 
lent width ratios toward p Oph A and C, Oph would result 
in 6-values of 0.60 km s^^ and 0.88 km s^^, respectively. 
Even if the uncertainties give infinite upper bounds, the 
nominal results are similar to the results of many earlier 
investigations. For instance, Palazzi et al. (1992) find 
h = 0.84 ±0.22 km in the direction of p Oph A, while 
Pan et al. (2004) Hst a value of 0.7 km s^^. For ( Oph, 
Crane et al. (1986) give b = 0.88 ± 0.02 km s'S and 
Roth & Meyer (1995) find h = 0.89 ± 0.19 km s'^. The 
CN absorption profiles toward HD 161056 arise from two 
components of approximately equal strength, separated 
by 5.9 km s^^. Interestingly, the i?(l)/_P(l) ratio for the 
component at —3.1 km s^^ indicates that this component 
has a slightly larger optical depth than the component 
at +2.8 km s-\ although the (0, 0)/(l, 0) R{Q) ratios in 
both cases exceed the optically thin limit. 

The UVES sight Hues in Table 5 other than HD 152236 
and HD 161056 exhibit fairly significant optical depths 
in their CN absorption lines. An extreme case is pre- 
sented by the sight line to HD 169454, for which the 
(0, 0)/(l, 0) i?(0) equivalent width ratio is only 2.83 ± 
0.05. Where measurements on both pairs of CN transi- 
tions arc available, the derived &- values agree remarkably 
well. Note that the ratio between the (0, 0) -R(O) and 
(1, 0) i?(0) equivalent widths provides the most strin- 
gent constraints on b because of the larger difference in 
intrinsic line strength for this pair, compared to J?(l) 
and P(l). The weighted mean of the two 6- values in 
the case of HD 154368 is 6 = 0.68 ± 0.01 km s'S while 
for HD 169454, HD 170740, and HD 210121 the weighted 
means arc b = 0.50±0.01 km s"\ 6 = 0.46±0.04 km s^\ 
and b = 0.80 ± 0.02 km s"'^, respectively. Palazzi et 
al. (1990) obtained a &-valuc of 0.65 ± 0.01 km s'^ for 
HD 154368, and the study by Roth & Meyer (1995) re- 
sulted in a 6- value of 0.68 ±0.08 km in this direction. 
Additionally, Palazzi et al. (1992) find b = 0.64 ± 0.17 
km s"^ for HD 169454, and b = 0.84 ± 0.20 km s'^ 
for HD 170740. Sheffer et al. (2008) give the 6-value 
toward HD 210121 as 0.9 km s~ . These various deter- 
minations show good agreement, except in the case of 
HD 170740. Yet, even for this sight line, the discrepancy 
is not that severe. The Palazzi et al. (1992) 6-valuc ac- 
tually agrees with our more precise determination within 
2 a and our result is particularly secure because it was 
derived from measurements of four different CN tran- 
sitions across two distinct bands, which were observed 
independently. Indeed, as a result of the excellent agree- 
ment between 6-values determined from the two equiv- 
alent width ratios, the weighted means for HD 154368, 
HD 169454, HD 170740, and HD 210121 were used as 
fixed input in the profile fitting analysis. For HD 73882, 
HD 152236, and HD 161056, the 6-values obtained from 
the R{1)/P{1) ratios were used as the initial input to 
ISMOD but were allowed to vary. 

The component parameters resulting from our profile 
synthesis fits to the CH+ and CN R{0) lines of the (0, 0) 
bands are presented in Table 6. In cases where isotopo- 
logic ratios were determined (see Table 7), the quoted 
velocities and 6-values apply to both the ^^C and ^•Re- 
hearing isotopologues since the two profiles were required 
to have identical component structure. Uncertainties in 
N for a given component were derived from the relative 
uncertainties in Wx. For 20 Tau and 23 Tau, the new 



results arc consistent with our previous analysis of CH"*" 
absorption in the Pleiades from McDonald spectra with 
much lower signal to noise (Ritchey et al. 2006). The 
only major difference is that, in the present investiga- 
tion, we derive a 10% higher column density of ^^CH+ 
toward 20 Tau, even though the equivalent widths agree 
to within 2%. As profile synthesis was not used pre- 
viously to obtain column densities, the discrepancy un- 
derscores the importance of carefully considering optical 
depth effects when dealing with strong absorption lines. 
For 23 Tau, both the equivalent widths and column den- 
sities of '^^CH"'" from this and the previous study arc con- 
sistent at the 1% level. Our new results on CH+ toward 
p Oph A are also broadly consistent with the analysis of 
Pan et al. (2004), who give a total ^^CH+ column den- 
sity that differs from our value by less than 1%. In the 
case of CN toward p Oph A, our profile fit results in a 6- 
value of 0.4 km s~^, somewhat lower than the 0.7 km s~^ 
given by Pan et al. (2004). Since the equivalent width 
measurements are essentially identical, the difference in 
b amounts to a 20% increase in A''(^^CN) in the present 
analysis. 

Table 7 lists the i^CN/^^CN and ^"^01+ /^^CB+ col- 
umn density ratios for our sample along with the iso- 
tope shifts and total molecular column densities de- 
rived through profile synthesis of the (0, 0) bands in 
CN and CH+. The uncertainties in i^CN/^^CN and 
^^CH+/^^CH+ were obtained by propagating the errors 
associated with the column densities of the two isotopo- 
logues, although the errors in A'' for the weaker isotopo- 
logue essentially determine the outcome. Also included 
in Table 7 are column densities and 12-to-13 ratios for 
CO from the literature, where such measm-cmcnts exist. 
In Table 8, we present a comparison between the col- 
umn densities that we derive from the (0, 0) R{0) and 
(1, 0) i?(0) fines of CH+ and CN. The good agreement 
indicates that we have adequately accounted for optical 
depth effects in these spectra, giving us confidence in our 
isotopologic results. 

As mentioned in § 2.2, we are not confident that 
^^CH+/^^CH+ ratios can be rehably obtained for the 
UVES sight lines. A number of contributing factors lead 
us to this assessment. There is a non- negligible amount 
of fringing in the reduced UVES spectra, making an un- 
ambiguous placement of the continuum impossible. The 
problem is compounded by the broad and moderately 
complex CH+ absorption profiles, which often consist 
of three or more components (see Table 6). Particu- 
larly worrisome are weak ^^CH"'' components that lie to 
the blue of the main component (for instance, toward 
HD 152236, HD 154368, and HD 169454). In the (0, 0) 
line at 4232 A, these components are likely to be blended 
with the strongest feature in ^■'CH+, considering that 
the isotope shift is —19 km s^^. In principle, one could 
constrain the ^^CH+ absorption profile using the (1, 0) 
transition at 3957 A since the ^^CH+ line corresponding 
to this transition is shifted by 33 km to the red of 
^^CH+, rather than to the blue. In practice, however, 
components that are already weak at 4232 A may not 
even be detected in the weaker transition at 3957 A, lim- 
iting the usefulness of this approach. 

Ultimately, wc did synthesize the CH"*" profiles toward 
the UVES targets in an attempt to derive ^'^CE.+ /"CH+ 
ratios, but the results bore out our initial concerns. No 
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Figure 12. Profile synthesis fits to tiie B—X (0, 0) band of CN toward HD 73882. See Figure 3 for a description of the plotting symbols. 
The same range in velocity is shown for each panel. The positions of the R{0) and R{1) lines of -'^CN are indicated by tick marks. When 
both isotopologues are present, the two profiles are synthesized simultaneously. 
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Figure 13. Same as Figure 12 except for the CN B—X (0, 0) band toward HD 152236. Two line-of-sight components are clearly identified 
in the R(0) and R(l) lines of ^•^CN and are also detected in the P(l) line. The position of the ^"^CN R{0) line corresponding to the stronger 
of the two components is indicated by a tick mark. Since the R{2) and P(2) lines are below the detection limit, no fits were attempted for 
these features. 



^^CH+ absorption is detected toward HD 152236 nor to- 
ward HD 210121 due to the low column densities of the 
individual ^^CH+ components in these directions. For 
HD 154368, HD 161056, and HD 169454, the ^'^CR+ ab- 
sorption profiles at 4232 A arc almost certainly blended 
with ^■'CH^, making it difficult to obtain any meaning- 
ful results on the ^^CH+/^^CH+ ratio. For sight lines 
with simpler profiles, the results are more encouraging, 
but the uncertainties are still quite large. The ^^CH"*" 
and 13CH+ features in the spectrum of HD 170740 do 
not appear to be blended, but ^■^CH+ is only detected 
at the 1.4-a level. The sight line to HD 73882 permits 
a relatively clean detection of both ^^CH+ and ^^CH+ 
at 4232 A, though the continuum level is somewhat un- 
certain. Nevertheless, we find a ^^CH+/^^CH+ ratio in 
this direction of 61.7 ± 17.3 for a fixed isotope shift of 
Av = -18.9 km s^^ 

We also analyzed the available UVES data on CH+ 
toward 23 Tau and ( Oph, both of which have simple, 
well-studied CH+ absorption profiles. From the (0, 0) 
transition at 4232 A, we obtained a ^'^CIl+ /^^CR+ ratio 
of 102.3 ± 38.4 for 23 Tau (Av = -19.4 km s'^) and 
88.5 ± 11.4 for ( Oph {Av = -19.0 km s'^). However, 
independent fits to the (1,0) line at 3957 A yielded lower 
values (55.5 ± 29.4 for 23 Tau and 52.8 ± 9.1 for C Oph). 



Taking the weighted mean of the two results in each case 
gives 72.8 ± 23.3 for 23 Tau and 66.7 ± 7.1 for ( Oph, 
both of which are consistent with the ^^CH+/^^CH+ ra- 
tios we derive from our McDonald observations of these 
stars (see Table 7). Still, the disparity in the values ob- 
tained from the two transitions liighlights the problems 
encountered when attempting to fit the undulating con- 
tinuum in the UVES spectra examined here. Given the 
uncertainties, our UVES results are not altogether dis- 
similar to those of Casassus et al. (2005) and Stahl et 
al. (2008), who first analyzed these data. Casassus et al. 
found a 12CH+/"CH+ ratio of 78.5±3.5 for the sight line 
to C Oph, while Stahl et al., reanalyzing the ( Oph data, 
obtained a ratio of 80.9 ± 3.0. In the direction of 23 Tau, 
Stahl et al. found i2CH+/i3CH+ = 96.5 ±16.3, while to- 
ward HD 73882, these authors obtained a i2CH+/i3CH+ 
ratio of 79.7 ± 14.5. 

The weaker lines of the B—X (0, 0) band in CN were 
fitted using the component structure found for the R{Q) 
line. Since the majority of CN absorption profiles consist 
of single velocity components, this restriction essentially 
amounted to keeping the 6- values fixed. In cases where an 
additional component is observed, the velocity difference 
between the two components, as determined from the 
R{0) line, was also held constant. The strengths of the 
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Figure 14. Same as Figure 12 except for tlie CN B-X (0, 0) band toward HD 154368. The positions of the R{0) and R{1) hnes of ^^CN 
are indicated by tick marks. The weak absorption feature blucward of the strong ^^CN R(0) Une is due to a second very weaJs component 
along the line of sight, which is positively identified in the spectra of CH+ and CH. 
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Figure 15. Same as Figure 12 except for the CN B—X (0, 0) band toward HD 161056. Two components of approxiiiiatcly equal strength 
are identified in each of the ^"^CN lines. The positions of the ^"^CN R{0) lines corresponding to these components are indicated by tick 
marks. 



individual components were allovifcd to vary, however, so 
that excitation temperatures could be derived separately 
for each cloud along the line of sight. Although two lines 
arising from the same excitation state, such as R{1) and 
P(l) or R{2) and P{2), should yield identical column 
densities, each line was synthesized independently in our 
analysis. The final column density for a given = 1 or 
N = 2 level was then determined by taking the weighted 
mean of the two results. Table 9 presents the CN column 
densities for each rotational level with a detectable tran- 
sition. For the four UVES sight lines where -^^CN i?(l) is 
detected, the R{1) lines of both isotopologues were fitted 
simultaneously as in the case of i?,(0). The same tech- 
nique was applied to the P(l) lines of ^^CN and ^^CN 
in the direction of HD 169454. Because this procedure 
permits the 12-to-13 ratios in different N levels to vary, 
we can accomodate small variations in excitation tem- 
perature between ^^CN and ^"''CN. Figiues 12 through 
18 give the profile synthesis fits to the CN B~X (0, 0) 
band for each of the UVES sight lines. 

Since the isotope shifts in CN and CH+ were left as 
free parameters during profile synthesis, it is instruc- 
tive to compare our results, from astronomical spectra, 
with those obtained in laboratory experiments. The 
mean isotope shift in our synthesis fits of the CII+ 



(0, 0) transition at 4232 A (toward McDonald targets) 
is At; = -18.9 ± 0.4 km s"! {5\ = -0.267 ± 0.005 A). 
This result is consistent with the experimental value of 
-19.3 ± 0.3 km s"! (-0.272 ± 0.004 A), which is based 
on the laboratory measurements of Carrington & Ram- 
say (1982) for and Benibenek (1997) for ^^GB+ . 
Similarly, our mean isotope shift for the CN (0, 0) i?(0) 
line, Aw = +13.1 ± 0.4 km s'^ {5\ = +0.169 ± 0.005 
A), agrees favorably with the laboratory value, Av = 
+13.6 km s"^ {5\ = +0.176 A; Bakker & Lambert 
1998). The four detections of "CN R{1) yield a mean 
isotope shift for the R{1) hne of +14.6 ± 0.4 km s~^ 
(+0.188 + 0.005 A), while, from the ^^CN and ^'''CN P(l) 
lines observed toward HD 169454, we derive an isotope 
shift of +9.2 ± 0.4 km s'^ (+0.118 ± 0.005 A). The lab- 
oratory results give +14.8 km s~^ (+0.192 A) for the 
R{1) line blend and +10.0 km s'^ (+0.130 A) for the 
P(l) blend (Bakker & Lambert 1998). The generally 
good agreement between the laboratory and astronomi- 
cal results serves as confirmation that our detections of 
the i^CN i?(l) and P(l) lines are real. 
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Figure 17. Same as Figure 12 except for tlie CN B-X (0, 0) band toward HD 170740. The position of the ^^CN R(0) line is indicated 
by a tick mark. 



4.1. ^'^CH+/^^CH+ Ratios 

One of the primary motivations for this study was 
to evaluate whether the reported variations in the 
^^CH+/^^CH+ ratio toward stars in the solar neighbor- 
hood are indicative of intrinsic scatter in the local inter- 
stellar ^^C/^^C ratio or result instead from systematic 
effects in the individual investigations. The results of 
our detailed analysis suggest that the latter scenario is 
more likely. In particular, we do not confirm the low ra- 
tios that were obtained by Hawkins & Jura (1987) from 
Lick Observatory data toward the Pleiades stars 20 Tau 
and 23 Tau (40 ± 9 and 41 ± 9, respectively), nor an 
earlier result toward 20 Tau (49^g^) by Vanden Bout & 
Snell (1980). Hawkins et al. (1985, 1989) reported low 
^2CH+/^^CH+ ratios based on Lick data toward C Oph, 
as well, giving values of 43 ±6 (Hawkins et al. 1985) and 
49± 10 (Hawkins et al. 1989). However, many subsequent 
investigations revealed ^^CH+/^'*CH+ ratios closer to 70 
for this line of sight, including Hawkins et al. (1993), 
who examined Kitt Peak National Observatory (KPNO) 
spectra of C Oph, arriving at a ratio of 63 ± 8. More 
precise results are available from Crane et al. (1991), 
who used the 2.7 m telescope at McDonald Observatory 
to derive a ^"^011+ /^^CB+ ratio of 67.6 ± 4.5, and Stahl 
& Wilson (1992), who reanalyzed ESO observations of 
C Oph (Stahl et al. 1989) and found a value of 71 ± 3. 



The Crane et al. (1991) result is, in fact, very close to 
the mean value of all of the C Oph results reported in the 
literature. Our new determination from McDonald data, 
while not as precise as the Crane et al. (1991) value, is 
nonetheless consistent with it. 

Hawkins et al. (1993) believe that the discrepancies 
between their KPNO results and the earlier Lick results 
for ( Oph (e.g., Hawkins et al. 1985) were caused by 
telluric contamination of the Lick spectra, which com- 
plicated the placement of the continuum. While this ex- 
planation may have resolved the controversy surrounding 
the 12CH+/13CH+ ratio toward C Oph, it left open the 
possibility that the spectra obtained by Hawkins & Jura 
(1987), including for 20 Tau and 23 Tau, were not sim- 
ilarly contaminated and that the low ^^CH+/^''CH+ ra- 
tios in these directions would later be confirmed. With 
our new results on ^2CH+/^^CH+ toward 20 Tau and 
23 Tau, such a position would now be difficult to sup- 
port. The Stahl et al. (2008) value for i2cH+/"CH+ 
in the direction of 23 Tau (96.5 ± 16.3) is further con- 
firmation that a ratio closer to 40 in the Pleiades clouds 
is unlikely, although the UVES results have their own 
complications (see § 3). Initially, we suspected that the 
difference between our McDonald results and those of 
Hawkins & Jura (1987) for the Pleiades stars was re- 
lated to the smaller 6-values we derive when modeling 
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Figure 18. Same as Figure 12 except for tlie CN B-X (0, 0) band toward HD 210121. The positions of the R(&) and R{1) lines of ^^CN 
are indicated by tick marlfs. 



the CH+ profiles. However, our optical depth correc- 
tions for these sight lines, as determined by the profile 
fitting routine, would only amount to increases of about 
4% in 7V(^^CH+) over the column densities obtained by 
Hawkins & Jura (1987) from their adopted corrections, 
whereas the ^^CH+/^^CH+ ratios differ by a factor of 
2. The dominant cause of the discrepancy is the almost 
factor-of-2 larger equivalent widths that Hawkins & Jura 
(1987) find for the ^^CB+ hues (see Table 2), which again 
points to possible errors in continuum normalization or 
to the effects of blending in their more moderate resolu- 
tion spectra. 

Federman et al. (2003) analyzed McDonald 2.7 m data 
toward p Oph A and found a ^"^CYi^ /^^GB+ ratio of 
120 ± 54, which, despite the large value, is within 1 cr of 
our more precise determination. These authors achieved 
a S/N ratio per pixel of about 650 in the vicinity of CH+ 
A4232, whereas the S/N is about 950 per pixel in our 
p Oph A data. As our spectra were also sampled at 
slightly higher resolution, this case illustrates the gen- 
eral trend that higher-resolution, higher signal-to-noise 
observations tend to result in ^^CH+/^^CH+ ratios that 
are closer to the average interstellar value of ^^G/^^C 
Indeed, our McDonald observations of 20 Tau, 23 Tau, 
p Oph A, and C Oph yield ^'^CB+ /^^CB+ ratios that 
are indistinguishable from 70 in each case. The weighted 
mean of the four determinations is 74.4 ± 7.6, or only 
0.6 a higher than the ^^C/^^C ratio adopted for the lo- 
cal interstellar medium. While this mean represents only 
a small sample of interstellar sight lines, the result is 
significant in the sense that each of these sight lines at 
one time had been found to have a much lower or much 
higher ^^CH+/^'^CH+ ratio. The statistical dispersion 
(i.e., standard deviation) in our new determinations of 
^2CH+/^3CH+ for these four sight lines (7.3) is essen- 
tially identical to the error in the weighted mean (7.6). 
We interpret these findings as confirmation of the theo- 
retical expectation that CH+/'^'^CH+ reflects the am- 
bient carbon isotopic ratio, and that this ratio does not 
vary within the solar neighborhood. 

4.2. ^^CN/^^CN Ratios 

The i^CN/i-'^CN ratios presented in Table 7 constitute 
the most extensive isotopologic results to date for CN in 
diffuse molecular clouds. The UVES results, in partic- 



ular, are considerably more precise than any literature 

results currently available. Gredel et al. (1991) reported 
abundance ratios of ^^CN to ^^CN for four interstellar 
sight lines. These included HD 73882 and HD 169454, 
for which they found i^CN/^^CN ratios of 170 ± 120 
and 38 ± 16, respectively. Due to the large uncertain- 
ties, these authors were unable to conclude whether or 
not the scatter in their sample represented true cloud-to- 
cloud variations. Our new results for these lines of sight 
reveal ^^CN/-'^'^CN ratios closer to the ambient value of 
i2Qyi3(^ |-,q|;1i instances. However, the small uncer- 
tainties in our determinations allow even modest varia- 
tions in the relative abundances of ^^CN and ^^CN to be 
discerned. Thus, the moderate increase in -'^CN/^'^CN 
that we find toward HD 73882 likely indicates an actual 
deficit in ^^CN, relative to ^^CN, for the gas in this direc- 
tion. Further confirmation of our result comes from an 
independent assessment by D. Welty (private communi- 
cation), who finds a i^CN/^^cN ratio of 87 for the line of 
sight to HD 73882. For HD 154368, Palazzi et al. (1990) 
reported a ^^CN/^^CN ratio of 101 ± 12, whereas we find 
a value essentially equal to 70. In this case, and in the 
comparison with the Gredel et al. (1991) results, the 
discrepancy in the ^^CN/^^CN ratio is largely driven by 
differences in the measured equivalent width of the ^^CN 
line (see Table 3). However, the Gredel et al. (1991) val- 
ues of Wa(^^CN) only differ from ours at the 1-cr level, be- 
cause of their significantly larger uncertainties, whereas 
the Palazzi et al. (1990) measurement of Wa(^^CN) for 
HD 154368 (0.82 ± 0.10 mA) is lower than our more pre- 
cise result (1.13 ± 0.04 mA) by 3 a. Equivalent width 
measurements of very weak lines are particularly sensi- 
tive to the placement of the continuum. However, by ex- 
amining the raw (unnormalized) spectrum of HD 154368, 
we can find no acceptable continuum fit that would yield 
a value of Wa(^^CN) as low as the Palazzi et al. (1990) 
result. 

Of the McDonald sight lines, both C Oph and C Per 
have previously been investigated for ^^CN/^'^CN. The 
most precise result available is that of Crane & Hegyi 
(1988), who found a i^CN/^CN ratio of 47.3^^^ in the 
direction of C Oph. Later, Roth & Meyer (1995) obtained 
a ratio of 35±13 for this line of sight. Our new determina- 
tion (48.8 ±19.5) has a precision comparable to the Roth 
& Meyer (1995) value, yet agrees well with the result of 
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Crane & Hegyi (1988). For C Per Kaiser et al. (1991) 
found a ^^CN/^^CN ratio of 77tll, wliicli is consistent 
with the shghtly less precise value we derive from our Mc- 
Donald observations of this star (67.0 ± 28.0). In order 
to obtain a more precise result for ^ Per, we can combine 
our determination with that of Kaiser et al. (1991) to 
arrive at a weighted mean ^^CN/^^CN ratio of 72 ± 19. 
Thus, again, we find a value of ^^CN/^'^CN that is in- 
distinguishable from the ambient ^^C/^'^C ratio. This 
is in contrast to the collection of ( Oph results, which 
provides strong evidence for a ^^CN/^^CN ratio that is 
significantly lower than the ambient value. We have al- 
ready noted the elevated ^^CN/^-'^CN ratio we find in the 
direction of HD 73882, which is significant because of the 
small uncertainty. Meyer et al. (1989) find a i^CN/^^CN 
ratio of 122 ± 33 toward HD 21483, which, despite the 
larger uncertainty, also suggests a higher-than-ambient 
value. Together, these results indicate that, unlike CH+, 
CN is affected by fractionation in at least some diffuse 
environments. 

The weighted mean of all of our determinations of 
12CN/13CN toward McDonald and UVES targets is 
67.5 ± 1.0, which agrees well with the adopted value for 
the local interstellar ^"^C/^^C ratio. However, there is 
significant scatter in ^^CN/^^CN from one sight line to 
another and the relatively small uncertainties in our de- 
terminations suggest that the variation represents true 
differences among the environments being probed. From 
our sample of ^^CN/^^CN ratios, we find a dispersion of 
24.9, which, in contrast to the situation for CH+, is much 
greater than the error in the weighted mean (1.0). These 
results for CN are similar, statistically, to published re- 
sults for CO. For example, the weighted mean value of 
i^CO/^CO from the ratios tabulated by Sheffer et al. 
(2007) is 69.9 ± 1.9, while the dispersion in those mea- 
surements is 29.5. Again, this suggests a level of fraction- 
ation in CN (and CO) not observed in CH+. However, 
some caution on this point may be warranted. First, we 
are unable to determine reliable ^^CH+/^^CH+ ratios for 
the majority of the UVES sight lines, meaning that our 
sample of CH"*" ratios is considerably smaller than that 
for CN. Second, the CN ratios toward HD 161056 and 
HD 170740, the two sight lines contributing the most to 
the scatter in ^^CN/^CN, are somewhat uncertain. The 
^^CN and ^^CN features are blended in the spectrum of 
HD 161056 because of the velocity separation between 
the two line-of-sight components. This may also con- 
tribute to the somewhat higher isotope shift we find in 
this direction. For HD 170740, there is some uncertainty 
in the placement of the continuum due to what may be 
either telluric absorption or the result of instrumental 
fringing just to the red of the ^^CN feature, which is also 
quite weak. Without these two extreme results, the scat- 
ter in i^CN/i'^CN is reduced to 10.3, though this might 
still be significant given the small error in the weighted 
mean value. 

4.3. 12 CN versus CO/^^ CO 

Many of the targets selected for this investigation 
were chosen because they exhibit fractionated values 
of 12CO/13CO (sec Table 7). That is to say their 
""^^CO/i^CO ratios arc significantly above or below the 
ambient value of ^^C/^^C. Since CN and CO are pre- 
sumed to reside in the same portion of a diffuse molecular 
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cloud (Pan et al. 2005), the expectation is that CN will 
be fractionated in the opposite sense compared to CO. 
This is due to the fact that CO is the most abundant 
C-bearing molecule in the ISM (by 2—3 orders of magni- 
tude over CN), and should regulate the availability of ^C 
atoms for incorporation into other C-bearing molecules 
that are cospatial and must draw from the same carbon 
reservoir. Figure 19 plots i^CN/^^CN versus i^CO/^^CO 
for the seven sight lines where both ratios have been de- 
termined. The i^CN/i'^CN ratios are from this investiga- 
tion (Table 7) except in the case of Q Oph, where we have 
used the more precise determination by Crane & Hegyi 
(1988), and in the case of C, Per, where we have plot- 
ted the weighted mean of our result and that of Kaiser 
et al. (1991). The i^co/^CO ratios are either from 
our archival study using HST (Sheffer et al. 2007) or 
from the literature (Hanson et al. 1992; Sonnentrucker 
et al. 2007). Also indicated in the figure is the value 
we adopt for the local interstellar ^^C/^'^C ratio (70 ± 7; 
see Sheffer ct al. 2007), which divides the parameter 
space into four quadrants. Two of these quadrants are 
of particular interest: the upper left quadrant, where 
low values of ^^CO/^^CO are associated with high val- 
ues of -^^CNZ-'^^CN, and the lower right quadrant, where 
high 12C0/"C0 and low i^CN/^-'^CN ratios coexist. In 
molecular gas characterized by a low ^^CO/^'^CO ratio, 
isotopic charge exchange has presumably converted some 
of the i^CO into ^^CO, while also exchanging ^^0+ ions 
for i^C"*". Since the formation of CN is initiated by re- 
actions involving 0+ (e.g., Federman et al. 1994). the 
gas in an environment where CO is fractionated by ICE 
should be enhanced in ^^CN relative to ^'^CN. In con- 
trast, a high i^CO/i'^CO ratio indicates that the molecu- 
lar gas has been subjected to selective photodissociation, 
which preferentially destroys ^^CO when only ^^CO is 
self shielded. In regions where SPD is occurring, the dis- 
sociated "'^■^C atoms are quickly ionized and incorporated 
into the pathways that lead to CN, which should result 
in enhanced I'^CN production. 

The data plotted in Figure 19 are suggestive of the in- 
verse relationship between ^^CN/^^CN and ^^CO/^^CO 
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that is expected due to photochemical fractionation and 
the physical association of CN and CO in diffuse cloud 
cores. Virtually no data points populate the upper right 
or lower left quadrants of the figure where the fraction- 
ation of CN and CO occurs in the same directions. The 
two data points that do (just barely) lie in these regions 
arc for the sight lines to 20 Aql (in the lower left) and 
^ Per (in the upper right). These sight lines have large 
uncertainties associated with their ^^CN/^^CN ratios, 
which in both cases are consistent with the unfraction- 
ated value. Although the downward trend in ^^CN/^^CN 
with increasing "'^^CO/^'^^CO is apparent, the correlation 
between these two parameters is rather shallow. In addi- 
tion, the appearance of a trend is somewhat contingent 
upon the data points at the extremes. The sight line 
with the lowest ^^CO/^^CO ratio (and highest value of 
12CN/13CN) is HD 73882, while the highest values of 
i^CO/i^cQ (^^^^ lo^gg^ i^cjN^/iacjN^ j.jjtiog) fo^nd in 

the directions of p Oph A and ( Oph (in order of increas- 
ing i^cO/i^CO). All of the i^CN/^^CN ratios between 
these extremes are consistent with the ambient value. 
Yet the uncertainties in i^CN/^^CN toward HD 73882 
and C Oph are among the smallest in the sample, giving 
credence to the perceived trend. 

The severity of CN fractionation depends not only on 
the effects that alter the relative abundances of ^^CO 
and ^'^CO but also on the balance between CO and C+, 
the other major reservoir of carbon in diffuse gas. If CO 
and C"*" have roughly equivalent abundances, then any 
alteration in ^^CO/^'^CO will have an equal, but oppo- 
site, effect on ^^C+ /^^C+ and hence on i^CN/^^CN. If 
0+ dominates over CO, however, fractionation in CO 
will have less of an impact on ^^CN/^^CN because there 
will be a larger pool of 0+ ions that will maintain the 
ambient 12-to-13 ratio. The reverse is also true. When 
CO dominates, CO fractionation will have a greater in- 
fluence on the ratio of ^^CN to ^^CN. A detailed model- 
ing effort to track how the relative abundance of each 
isotope (and isotopologue) of interest depends on the 
physical conditions in the cloud and the parameters of 
the external UV radiation field is beyond the scope 
of this paper. Instead, we employ a simple, analyti- 
cal method to estimate the fraction of carbon that is 
locked up in CO along a given line of sight based on 
the measured isotopologic ratios in CN and CO. In what 
follows, we define /(^^CN) = iV(i-^CN)/[Ar(i2CN) + 
Ar(i3cN)] = l/[Af(i2CN)/iV(i3CN) + 1] as the fraction 
of CN molecules that take the form of ^''CN. Similar 
definitions can be made for ^'^CO and ^'^CH+, where 
/(13CH+) = /(i-'C) = 1/(70 + 1) - 0.014 is taken 
to be the ambient fraction of ^'^C. With these defini- 
tions, we can write /(^^CN) as a function of /(^^CO) 
and /("CH+). We have 

/(i3cN) = /(i3CH+)[l-x] 

+ {/(^^CH+) + [fC^CR+) - fC^CO)] } X 
= /(^3CH+) + [/(13CH+) - fC^CO)] X, (1) 

where the factor x determines how strongly /(-'^^CN) 
is affected by changes in f(}^CO) and is equivalent to 
the column density ratio iV(CO)/A^(C+), assuming the 
abundance of neutral carbon is negligible. In the limiting 
case a; = [i.e., Ar(C+) > Ar(CO)], CN is unaffected by 



fractionation in CO and ^^CN/^'^CN equals the ambient 
ratio. 

Equation (1) can be rearranged to yield an expression 
for X, 

^ /(^3CN)-/(13CH+) 

/(13CH+) -/(13C0)' ^ ' 

which is valid whether CO is fractionated upward via 
selective photodissociation or downward due to isotopic 
charge exchange. However, when selective photodisso- 
ciation is under consideration, the fractions should be 
defined as, for example, /(^^CN) = N(^^Qn)/N{^^Gn), 
because SPD docs not exchange one isotopologue for an- 
other but simply destroys the less abundant isotopologue 
leaving the more abundant one unchanged. The frac- 
tion of total carbon that is locked up in CO can be 
found from x by /(CO) = 7V(C0)/[iV(C0) + N{G+)] = 
[N{CO)/N{C+)]/[N{CO)/N{G+) + 1] = x/{x + 1). We 
use this relation in conjunction with equation (2) to de- 
rive /(CO) for the sight lines included in Figure 19. The 
resulting values are /(CO) = 0.09 for HD 73882, 0.01 
for HD 154368, 0.07 for HD 210121, 0.40 for p Oph A, 
and 0.45 for Q Oph. No physically meaningful determi- 
nations of /(CO) are possible toward 20 Aql or ^ Per 
because the measured isotopologic ratios in CN and CO 
fall on the same side of the ambient ratio. Again, we note 
that in these cases the ^^CN/^^CN ratios are consistent 
with being unfractionated. 

For the sight lines in Ophiuchus, we find that nearly 
50% of the carbon is contained in CO, which is consis- 
tent with the fact that the UV optical depths (tuv = 
1.98 and 2.04 for the clouds toward C, Oph and p Oph A, 
respectively) are near that where the C+ to CO con- 
version is expected to occur (i.e., ruv = 2; see Fed- 
crman et al. 1994). The value of /(CO) is substan- 
tially lower toward HD 210121 even though the UV op- 
tical depth in this direction (ruv = 3.35; Sheffer et 
al. 2008) is larger than those for the Ophiuchus clouds. 
Here, we note that the ^^CO/^^CO ratio in the direc- 
tion of HD 210121 is consistent with 70 and so does 
not show strong evidence of fractionation, meaning that 
/(CO) cannot be determined. An mifractionatcd value 
of i^CO/i^cO toward HD 210121 is also consistent with 
the suggestion by Welty & Fowler (1992) that the UV ra- 
diation field incident on the cloud in this direction (the 
high-latitude molecular cloud DBB 80; see de Vries & 
van Dishoeck 1988) is lower than the average interstel- 
lar field. Finally, the low values of /(CO) we find to- 
ward HD 73882 and HD 154368 are surprising consider- 
ing the much larger UV optical depths in these directions 
(adopting ruv — ~ 4.7 in both cases). The larger 
values of tuv indicate that the sight lines to HD 73882 
and HD 154368 probe deeper into their respective clouds, 
where a majority of the carbon is expected to be in the 
form of CO rather than in neutral atomic or ionic form. 
However, equation (2) only applies if the ^^CN/-^^CN ra- 
tio reflects the ^^C+/^^C+ ratio of the gas from which 
CN formed. If ^^CN/i^qjn^ substantially altered after 
CN formation, then equation (2) no longer holds. Thus, 
instead of surmising that there is a very low fraction of 
elemental carbon in CO in the lines of sight to HD 73882 
and HD 154368, we suggest that we are seeing the effects 
of a separate fractionation process in CN, which acts to 
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lower ^^CN/^"^CN such that the enhancements due to CO 
fractionation arc not as large as expected. 

We assume that the low ^^CO/-^'^CO ratios in the di- 
rections of HD 73882 and HD 154368 result from an ef- 
ficient isotopic charge exchange reaction involving CO, 
which implies a colder kinetic temperature for the molec- 
ular gas along these lines of sight. Yet, at colder tem- 
peratures, CN is similarly aifected by ICE via the reac- 
tion + 12CN ^'^C+ + i^CN + AE, where the 
zero-point energy difference AE/ks = 31 K (Bakker & 
Lambert 1998; derived from the molecular parameters in 
Prasad ct al. 1992) favors ^^CN. Unlike CO, however, 
CN is dissociated by continuum radiation, not by lines, 
and so isotope-selective self shielding should not be a fac- 
tor. Any dissociating photons will equally destroy ^^CN 
as well as ^^CN. Yet if the timescale for the ICE reaction 
with CN is shorter than the photodissociation timescale, 
then the reaction will reach equilibrium and we can write 



exp 



AE 



n(i^CN) _ n(i^C+) kr _ nC^C+) 
n(i3CN) ~ n(i3C+) kf ~ n(^^C+) 

where n{X) represents the number density of species X, 
kf is the rate coefficient of the forward reaction, and kr = 
A;/exp(— Ai^/fcBTkin) is the rate coefficient of the reverse 
reaction. Analogous expressions have been derived in the 
case of CO (see Sheffer et al. 1992; Lambert ct al. 1994). 
When photodissociation of CN dominates over isotopic 
charge exchange, no upward fractionation occurs, as it 
does for CO, because the photodissociation rates of the 
two CN isotopologues are equal (i.e., = Tis) and we 
have 



n(i^CN) _ ri(i2c+) £13 _ ni^^C+) 
n(i3CN) ~ n(i3C+)ri2 " n{^^C+)' 



(4) 



It is important to note that the ^'^C'^ /^^C~^ ratio in 
equations (3) and (4) refers to the ratio that has been 
elevated above the ambient value of ^"^C/^^C by isotopic 
charge exchange with CO. Thus, when isotopic charge 
exchange with CN lowers ^^CN/^'^cn, this ratio is still 
higher than the ambient value. If we adopt a kinetic 
temperature of 51 K for the clouds toward HD 73882 
and HD 154368, a value which corresponds to the ob- 
served H2 rotational temperature (Tqi) for both sight 
lines (Rachford et al. 2002), then we find, from equa- 
tion (3), that ICE would have lowered i^CN/^^CN by 
a factor of 1.84. If the actual ^^C^ /^^C^ ratios toward 
HD 73882 and HD 154368 are higher than the measured 
^^CN/^^CN ratios by this amount, then the values of 
/(CO) in these directions would be 0.24 and 0.34, re- 
spectively. However, the rotational excitation tempera- 
ture in C2 (71)2) may provide a better indication of the 
kinetic temperature in the denser central portion of a 
diffuse cloud where CN resides. Indeed, Sonnentrucker 
et al. (2007) find that To2(C2) is consistently lower than 
2oi(H2) for a given sight line. These authors fit the- 
oretical models for the excitation to the observed C2 
rotational level populations to estimate the actual ki- 
netic temperature, finding Tkin — 20 ± 5 K for both 
HD 73882 and HD 154368. At this kinetic tempera- 
ture, the true ^^C+/^3c+ ratios would be a factor of 
4.71 higher than i^CN/^^CN, resulting in /(CO) values 
of 0.32 for HD 73882 and 0.48 for HD 154368. Thus, 



the inclusion of an isotopic charge exchange reaction 
with CN leads to more realistic values for the fraction 
of C in CO for sight lines with low ^^CO/^^cq ratios. 
However, even at 20 K, the derived value of /(CO) to- 
ward HD 73882 is smaller than would be expected. With 
our simple analytical model, it is difficult to accommo- 
date a -'^^CO/'^'^co ratio as low as 25 and also demand 
that a high fraction of carbon be locked up in CO. The 
uncertainty in the determination of ^^CO/^ 'CO toward 
HD 73882 is considerable, however, and indicates that a 
ratio as high as 47 would be consistent with the observa- 
tions. At Tkin = 20 K, a ^^CO/^CO ratio of 47 would 
result in a value of /(CO) = 0.63 for HD 73882. In the 
case of HD 154368, the maximum ^^co/^^co ratio al- 
lowed by the much smaller reported uncertainty is 45, 
which, at Tkin = 20 K, would yield /(CO) = 0.59. As 
these calculations demonstrate, there is a real need for 
expanding on existing chemical models of diffuse clouds 
in order to clarify the situation regarding photochemical 
fractionation in all of the relevant C-bearing molecules 
and atomic carbon species. 

5. CN ROTATIONAL EXCITATION 

The CN column densities resulting from our profile 
synthesis fits to the B—X (0, 0) band can be used to de- 
rive the CN rotational excitation temperatures Tqi and 
T12 by applying the Boltzmann equation 



hp 

fee 



In 



9j 



(5) 



where g and N denote the statistical weight and col- 
umn density of the lower (ith) or upper (jth) rota- 
tional level. The energy separation associated with the 
N = 1—0 transition at A = 2.64 mm is taken to be 
hv/k^ = 5.445 K, while for N = 2-1 at A = 1.32 mm, 
we assume hv/k-B = 10.885 K (Black & van Dishoeck 
1991). Uncertainties in T^- were calculated from the col- 
umn density uncertainties {(Jn) through standard error 
propagation such that 
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(6) 

The excitation temperatures and 1-cr errors obtained 

from equations (5) and (6) are listed in Table 9, along 
with the column densities of ^^CN in the first three rota- 
tional levels of the ground vibrational state. Also shown 
are the values of Tqi for the ^^CN molecule for the four 
UVES sight lines with detections of ^^CN R{1). For 
HD 73882, HD 154368, and HD 210121, the A^ = 1 col- 
umn density of -'^^CN is based solely on the value obtained 
from the i?(l) line. Since both the -R(l) and P(l) lines 
of i^CN are detected in the spectnun of HD 169454, the 
final value of N{N = 1) for ^'^CN in this direction is the 
weighted mean of the column densities derived from the 
two transitions. 

As expected, our derived rotational temperatures are 
consistent with the dominant excitation mechanism be- 
ing radiative pumping from the cosmic microwave back- 
ground (see Figure 20). The weighted mean value of 
Toi(^2CN) in our sample is 2.754±0.002 K, which implies 
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an excess over the temperature of the CMB (Tcmb = 
2.725 ± 0.002 K; Mather et aL 1999) of only 29 ± 3 mK. 
This modest excess hkely arises from local collisional 
excitation by electrons in some environments (see be- 
low). However, the dispersion in individual values of Tqi 
(134 mK) cautions that the excess may not be physical 
but rather due to measurement or modeling uncertain- 
ties in our analysis. There is a slight indication of a 
greater excess in Ti2(^^CN) based on the weighted mean 
of 2.847 ± 0.014 K, but the dispersion in these measure- 
ments (262 mK) is also greater. An elevated value of 
Ti2 could be caused by an overcstimation of the strength 
of the weak R{2) and P{2) lines, resulting from noise 
in the spectrum. This is likely to blame for the higher 
values of T12 we find toward the McDonald targets, as 
well as toward HD 161056 and HD 170740. These are 
the sight lines with the weakest CN absorption in our 
sample. For the directions where CN absorption is the 
strongest, namely, HD 73882, HD 154368, HD 169454, 
and HD 210121, the agreement between Tqi and T12 is 
at the 3% level or better. The rotational excitation tem- 
peratures derived from the observed ^^CN i?(0), -R(l), 
and P{1) lines show no excess over the temperature of 
the CMB. Quantitatively, the weighted mean value of 
Toi("CN) is 2.652 ± 0.069 K, which is within 1 a of 
Tcmb- 

Recently, a compilation of CN excitation temperatures 
along over 100 lines of sight was published by Slyk et 
al. (2008). These authors claim to detect an excess 
over the cosmic microwave background of 580 mK based 
on their average temperature of 3.31 K. Upon closer in- 
spection, however, the rather high values of Tqi given 
in Slyk et al. (2008) are almost certainly the result of 
not accounting adequately for optical depth effects in 
the strong i?(0) line. While the authors make numerous 
attempts at estimating the Doppler widths, they settle 
on a simple assumption that b = 1.0 km for all di- 
rections. In our much smaller sample, the vast majority 
of CN components have 6 < 1.0 km s~^ and we note 
that differences of only a few tenths of a km s~^ can 
greatly impact the derived column densities, even for 
lines with moderate optical depths. If we restrict the 
Slyk et al. (2008) sample to only the eleven sight lines 
in common with our investigation, then their average ex- 
citation temperature is still 3.30 K. The largest discrep- 
ancies between the two samples arise for the three sight 
lines with the most severe cases of optically thick absorp- 
tion, HD 73882, HD 154368. and HD 169454, for which 
Slyk et al. (2008) give Tm 3.68 K, 3.67 K, and 5.06 K, 
respectively. In each of these cases, their measured val- 
ues of Wx agree with ours at the 5% level. However, their 
column densities, especially those derived from the R{0) 
line, are significantly smaller than ours [by 40% for the 
i?(0) lines toward HD 73882 and HD 154368 and 76% in 
the case of i?(0) toward HD 169454]. The smaller values 
Slyk et al. infer for N{N = 0) are then directly respon- 
sible for their higher values of Tgi. 

Palazzi et al. (1992) previously examined CN rota- 
tional excitation in a sample of 34 interstellar sight lines, 
ten of which are included in our investigation. Their 
weighted mean of Tqi = 2.817 ± 0.022 K is within 3 a 
of our more precise result. From their mean value, 
Palazzi et al. (1992) deduced an excess temperature of 
82 ± 30 mK, which would be 92 mK relative to the value 
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Figure 20. Rotational excitation temperatures in ^•^CN versus 
tlie CN/CH+ eolumn density ratio for individual velocity compo- 
nents. The full 3-0" range in T()i(^^CN) is shown in each case. 
The solid line denotes the temperature of the Cosmic Microvirave 
Background (Tcmb = 2.725 K; Mather et al. 1999), while the 
dashed line represents the weighted mean value of Tqi for our sam- 
ple (2.754 K). The dotted lines indicate the (unweighted) mean and 
stajidard deviation for this set of measurements (2.707 ±0.134 K). 
The two data points that lie above the l-cr range for the sam- 
ple are for the main component toward HD 154368 and for the 
—3.1 km component toward HD 161056 (in order of increasing 
CN/CH+). 

of Tcmb adopted here. Again, if we consider only the 
ten sight lines from Palazzi et al. (1992) that are also 
analyzed in the present work, their average excitation 
temperature is 2.813 K, which is essentially identical to 
their full sample mean. Evidently, sample size has little 
impact on the mean value of Toi("'^^CN), which is more 
heavily dependent on systematic effects, and, specifically, 
on corrections for optical depth. The basic difference be- 
tween our study and that of Palazzi et al. (1992) is that, 
for the ten sight lines in common, Palazzi et al. find b- 
values that are, on average, 0.1 km s~^ higher than the 
ones we derive. This is a minor difference, but it ac- 
counts for the 2% higher mean excitation temperature 
that these authors obtained. We also note that, for this 
sample of ten sight lines, we have reduced the scatter in 
Toi from 0.3 K to 0.1 K. 

Although our mean CN excitation temperature 
(2.754 ± 0.002 K) implies a smaller excess over the tem- 
perature of the CMB (29 ± 3 mK) than has been found 
previously, our result has greater statistical significance 
because of the reduced uncertainties. However, two of 
the four sight lines in our sample with the strongest 
CN absorption lines, and therefore; the smallest errors 
in Tqi, have excitation temperatures below Tcmb- These 
are the sight lines to HD 73882, where Tqi = 2.631 ± 
0.004 K (3.4% below Tcmb), and HD 210121, for which 
Toi = 2.604 ± 0.007 K (4.4% below Tcmb)- The ro- 
tational excitation temperatures in the other two direc- 
tions with strong CN arc significantly above Tcmb- For 
HD 154368, we find Tqi = 2.911 ± 0.004 K (6.8% above 
Tcmb), while for the sight line to HD 169454, we obtain 
Toi = 2.804 ±0.003 K (2.9% above Tcmb)- Since there is 
no known physical mechanism that can lower Tqi below 
the temperature of the cosmic microwave background, 
our resuhs for HD 73882 and HD 210121 indicate that 
perhaps the uncertainties have been underestimated in 
these instances. The sight lines to HD 73882, HD 154368, 
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HD 169454, and HD 210121 exhibit the largest CN op- 
tical depths in our sample, meaning that the excitation 
temperatures are sensitive to small changes in b. In- 
creases in b of only 0.03 km and 0.06 km s~^ for 
HD 73882 and HD 210121, respectively, would be suffi- 
cient to bring the derived values of Tqi into agreement 
with TcMB- (Likewise, for HD 154368 and HD 169454, 
decreases in b of 0.04 km s~^ and 0.01 km s^^, respec- 
tively, would eliminate the small excesses in Tqi observed 
in these directions.) Such small changes would not be 
unreasonable considering the level of variation present in 
the 6-values obtained through the doublet ratio method 
for these sight lines (see Table 5). 

Figure 20 plots our derived CN excitation tempera- 
tures for individual velocity components (shown as 3-a 
error bars) against the CN/CH+ column density ratio, 
an empirical density indicator (see Sheffer et al. 2008). 
The CMB temperature (2.725 ± 0.002 K; Mather et al. 
1999) and our weighted mean value (2.754 ± 0.002 K) 
are indicated in the figure. Also shown are the (un- 
weighted) mean and standard deviation for our sample of 
excitation temperatures (2.707±0.134 K). In all but two 
cases, the derived values of Tqi are within one standard 
deviation of the unweighted mean, even if the individ- 
ual errors are much smaller. (The small error bars at 
high density that are below Tcmb arc for the sight lines 
to HD 73882 and HD 210121 discussed above.) Only 
in the main component toward HD 154368 and in the 
—3.1 km s~^ component toward HD 161056 is Tqi more 
than 134 mK larger than the average and both of these 
clouds have high CN/CH+ ratios. For the latter com- 
ponent, we find Tqi = 2.960 ± 0.018 K. This resuh may 
be particularly significant since the CN column densities 
toward HD 161056 seem not to be affected by optical 
depth effects in any appreciable way. Additionally, the 
excitation temperature in the +2.8 km component 
in this direction (Tqi = 2.761 ± 0.020 K) is consistent 
with Tcmb at the 2-a level and this component has a 
CN/CH+ ratio that is a factor 30 lower than the ratio in 
the —3.1 km s~^ component. A rotational temperature 
in excess of Tcmb results from a local excitation mecha- 
nism, presumably electron impact (Thaddeus 1972). As 
such, an excess in Tqi might be expected at higher gas 
densities, if the electron density in the gas is also higher. 

To further explore the local excitation process, we can 
use our determinations of Tqi to derive estimates for the 
electron density rig in these clouds. Our approach differs 
from previous analyses of this kind (e.g., Meyer & Jura 
1985; Palazzi et al'. 1992; Roth & Meyer 1995) in that 
past investigators have used an independent measure of 
the contribution to Tqi from local coUisional excitation, 
Tioc, to correct their Tqi values in order to derive the 
temperature of the CMB. These studies have estimated 
the local contribution either from ionization balance by 
using column densities of Ca I and Ca II, for example, 
to calculate rig, or from microwave observations of CN 
emission at 2.64 mm to evaluate Tioc directly. In our 
analysis, we adopt Tcmb = 2.725 K, and assume that 
any additional excitation results from electron collisions. 
The excess temperature, then, yields an estimate for the 
electron density in the gas. Following Thaddeus (1972), 
the rate of excitation from N = to N = 1 may be 
written 



where Ai = 1.19 x 10"^ s"^ (Thaddeus 1972) is the spon- 
taneous radiative decay rate of the iV = 1 level and Cqi 
represents the total rate, including both radiative and 
coUisional excitation. If Tqi = Tcmb, then the excitation 
rate is Coi = Ccmb = 5.60 x 10^^ s^^. Any local exci- 
tation will increase this rate. Thus, when Tqi > Tcmb, 
the local rate is given by Cioc = Cqi — Ccmb- Once 
Cioc is known, the electron density can be found from 
ne = Cioc/foi, where roi is the rate coefficient for CN 
rotational excitation by electron impact. 

We calculated rie from equation (7) and the above re- 
lations for clouds with significant excesses in Tqi by as- 
suming a gas kinetic temperature of 20 K. This temper- 
ature corresponds to that calculated by Sonnentrucker 
et al. (2007) for the clouds toward HD 154368 and 
HD 169454 from their analysis of C2 rotational exci- 
tation. The rate coefficient given by Allison & Dal- 
garno (1971) for CN excitation by electron impact at 
Tkin = 20 K is roi = 1.3 x 10"^ cm^ 3"^ For the main 
component toward HD 154368, which has an excess ro- 
tational temperature of 186 mK, these calculations yield 
an electron density of 0.69 cm~^. For the —3.1 km s~^ 
component toward HD 161056, where the measured ex- 
cess is 235 mK, we find rig = 0.88 cm^"^. Finally, for 
HD 169454, the excess of 79 mK gives = 0.29 cm-^. 
These electron densities are generally larger than those 
derived by Black & van Dishoeck (1991), who found val- 
ues of He in the range 0.02—0.5 cm^'^, and would indicate 
large gas densities in the cloud cores (uh ~ 6000 cm~^ 
for the component at —3.1 km s~^ toward HD 161056, 
assuming ne/nn ~ 1-4 x 10""^). The problem is even 
more severe in cases where /(CO) is substantial [e.g., 
toward HD 154368, where we find /(CO) 0.5], since 
fewer carbon atoms are available to provide free elec- 
trons. For comparison, Sonnentrucker et al. (2007) find 
riH ~ 250 cm-3 for HD 154368 and HD 169454 based 
on their C2 analysis, while Sheffer et al. (2008) ob- 
tain riH ^ 750 cm~^ for the main component toward 
HD 154368 using an analytical expression involving the 
formation and destruction pathways for CN. 

It is probably more appropriate to consider the electron 
densities derived from the excess rotational temperatures 
in CN to be upper limits because excitation from neu- 
tral impacts (mainly H2) can also contribute. Moreover, 
the actual excess temperature due to local coUisional ex- 
citation is likely to be less than the measured value of 
Tioc = Tqi — Tcmb, considering the dispersion of 134 mK 
in our determinations of Tqi. Taking this dispersion into 
account, the actual excess temperatures for the clouds 
toward HD 154368 and HD 161056 might be as low as 
52 mK and 101 mK, respectively. Indeed, from observa- 
tions of CN emission at 2.64 mm, Palazzi et al. (1990) 
find Tioc = 35 ± 9 mK for the main component toward 
HD 154368. These lower excesses would yield electron 
densities of rie = 0.19 cm'^ for HD 154368 and 0.37 cm'^ 
for the —3.1 km s~^ component toward HD 161056. An 
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excess in Toi would then not be inferred for HD 169454, 
since the observed rotational temperature in this direc- 
tion is within the typical dispersion of the sample. Here, 
again, we see a need for improved detailed models, in 
this case, of the excitation processes in diffuse molecular 
clouds. The precise set of observational results presented 
in this work would offer robust constraints for any such 
modeling efforts. 

6. CONCLUDING REMARKS 

Isotopologic ratios in interstellar molecules containing 

carbon are effective probes of the chemical processes ac- 
tive in diffuse environments. In this investigation, we 
examined optical absorption lines of CN and CH+ along 
a total of thirteen lines of sight through diffuse molec- 
ular clouds. The very high signal-to-noise ratio obser- 
vations enabled us to extract precise ^^CN/^^CN and 
^^CH+/^^CH+ ratios, which were used to assess various 
predictions of diffuse cloud chemistry. 

Our results on ^^CII+/^^CII+ confirm that this ratio 
does not deviate from the ambient ^^C/^'^C ratio in lo- 
cal interstellar clouds, as expected if CII+ is formed via 
nonthermal processes. Each of the McDonald sight lines 
in our sample with detectable absorption from ^"^011+ 
had at one time been found to have cither a much lower 
or a much higher value of ^^CH^/^'CH+. Such values 
are not confirmed here. In particular, the ^^CH+/^^CH+ 
ratios that we derive toward the Pleiades stars, 20 Tau 
and 23 Tau, are both indistinguishable from 70, preclud- 
ing a ^^C/^^C ratio near 40, as advocated by Hawkins 
& Jura (1987), for the gas associated with this nearby 
cluster. Considering these results, we find no evidence 
for variation in ^^C/^'^C within the solar neighborhood. 
In contrast, the studies by Casassus et al. (2005) and 
Stahl et al. (2008), based on UVES spectra, found vari- 
ation in 12CH+/13CH+ at about the 16% level. These 
authors attributed the variation to intrinsic scatter in 
the local ^^C/^^C ratio due to inefficient or incomplete 
mixing in the interstellar medium. From our analysis of 
some of the same UVES data, we conclude that the vari- 
ation in 12CH+/"CH+ is likely related to the difficulties 
of obtaining accurate ^^CH+/^^CH+ ratios from these 
spectra. The major difficulties involve the placement of 
the rather poorly-defined continua, and the blending of 
absorption from ^^CH"'" with weak ^^CH+ components 
blueward of the main components. It is probably not a 
coincidence that the ^^CH+/^^CH+ ratios derived from 
our high-resolution McDonald spectra, which have flat, 
well-behaved continua, along lines of sight with simple 
CH+ absorption profiles show little variation. While in- 
trinsic scatter in the local ^^C/^'^C ratio may exist at 
some level, due to small differences in the amount of 
chemical enrichment and subsequent mixing from one 
site to another, it is probably too early to say that this 
variation has been detected. 

Unlike in the case of CH+, the isotopologic ratios in 
CN and CO show evidence for significant fractionation 
away from the ambient value of ^^C/^^C in some dif- 
fuse environments. We attribute these results to the 
competing influences of isotope-selective photodissocia- 
tion and isotopic charge exchange reactions, which vary 
with the strength of the UV radiation fleld and the physi- 
cal conditions within a cloud. The results on ^^CN/^'^CN 
and ^^CO/^^CO in lines of sight where both ratios have 



been determined are suggestive of the inverse relation- 
ship between these quantities that is anticipated, assum- 
ing that CN and CO coexist in diffuse cloud cores. This 
relationship was flrst directly suggested by Sheffer et al. 
(2007), but the isotopologic results for CN that existed 
at the time were not extensive enough to adequately test 
the idea. With the results presented here, the situa- 
tion is now significantly improved, but further observa- 
tions of both i^CN/^CN and ^^CO/^CO would help 
to clarify the precise relationship between the two quan- 
tities. Specifically, measurements of ^^CO/^^CO toward 
HD 169454, where i^CN/^^CN is unfractionated, and to- 
ward HD 170740, where CN is fractionated upward (by 
a factor of 1.9 ±0.5 over the ambient value), would prove 
useful. One aspect of the chemistry of diffuse clouds that 
future modeling efforts will need to examine is the role 
played by isotopic charge exchange with CN, which we 
have found to be important in cold diffuse cloud cores. 
An ICE reaction involving CN is the most logical ex- 
planation for why the enhancements in ^^CN/^^CN are 
not as large as expected along lines of sight with low 
i2C0/"C0 ratios, such as HD 73882 and HD 154368. 

The effects of photochemical fractionation in carbon- 
bearing molecules are predicted to be the strongest in 
diffiise moleciilar clouds becaiise dark clouds are better 
shicldcid from UV radiation and lack a significant abun- 
dance of carbon ions. This has been convincingly demon- 
strated by millimeter-wave emission studies of CN, CO, 
and H2CO (see Milam et al. 2005). However, it is in- 
teresting to note that the expression developed in § 4.3, 
which yields the fraction of carbon locked up in CO, may 
still apply in denser regions. Keene et al. (1998) de- 
tected 809 CHz emission from both ^^C" and ^^C° near 
the center of the Orion Nebula, finding a i2qo^i3(^o 
tio of 58 ± 12. This was after Boreiko & Betz (1996) 
observed 158 /im emission from both isotopic variants of 
C+ in this region and found a ^^C^/^^C^ ratio of 58 ±6. 
The i2ci8o/i3ci80 ratio, also derived by Keene et al. 
(1998) from observations of the 2—1 and 3—2 transitions, 
is 75 ± 9. While all of these ratios are consistent with an 
ambient ^^C/^^C ratio of 70, the exact correspondence 
between the ratios in neutral and ionized carbon sug- 
gests that there is a real difference between those ratios 
and that found in C^*0. Since these observations sam- 
ple gas in a well-known photon-dominated region (PDR) , 
it is reasonable to assume that the slight increase in the 
12q18q^13qisq j-g^^j^ ^y^j. ^j^g ambient value is the result 

of selective photodissociation. An application of equa- 
tion (2) to these measurements, with the substitution of 
/(i'^C+) for /(i^CN), would yield a value of /(CO) = 
0.76, which is consistent with expectations for gas in a 
dense molecular cloud like that in Orion. We contrast 
this result with the values of /(CO) we obtain toward 
the diffuse molecular clouds in Ophiuchus, namely 0.40 
toward p Oph A and 0.45 toward ( Oph. The Ophiuchus 
results are appropriate for sight lines that probe the tran- 
sition region between diffuse and dense clouds, where the 
conversion of C+ into CO is expected to occur. 

Beyond isotopologic ratios, the high quality optical 
spectra examined in this investigation allowed us to ob- 
tain precise rotational excitation temperatures for both 
i^CN and i^CN. Our weighted mean value oiToi{^^CN) 
= 2.754±0.002 K implies an excess over the temperature 
of the cosmic microwave background of only 29 ± 3 mK. 
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Figure 21. Synthetic fits to the Ca II K profiles toward the unred- 
dened stars a Leo {upper panel) and a Vir {lower panel). Plotting 
symbols are the same as in Figure 3. The Ca II profile in the di- 
rection of a Leo is synthesized using two line-of-sight components, 
with LSR velocities of +4.4 and +9.7 km s^^. For a Vir, five 
components comprise the fit. These have LSR velocities of —12.7, 
-8.9, -6.4, -1.5, and +0.5 km s"!. 

This excess is considerably smaller than that suggested 
in the recent survey by Slyk et al. (2008), but is also re- 
duced compared to the excess found by Palazzi et al. 
(1992). Both of these previous studies indicated the 
need for an additional excitation mechanism beyond elec- 
tron and neutral collisions (and the CMB) to account 
for the observed excitation of CN. The more modest ex- 
cess found here eliminates the need for such a mecha- 
nism in the general ISM. Only three sight lines in our 
sample exhibit significant excesses in Tqi over the cos- 
mic microwave background temperature. For the clouds 
in these directions, we derived upper limits to the elec- 
tron densities in the range 0.3—0.9 cm"*^. The rotational 
excitation temperatures observed in ^^CN, from strong 
detections of the R{0) and R{1) lines in four directions 
and the first interstellar detection of the P(l) line, show 
no excess over the temperature of the CMB. 
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APPENDIX 

INTERSTELLAR ABSORPTION TOWARD 
UNREDDENED STARS 

In this appendix, we provide results on weak Ca I 



(A4226) and Ca II K (A3933) absorption toward the stars 
a Leo and a Vir, which served as unrcddcncd standard 
stars for our McDonald observations. Although a Lyr 
was also observed once as a standard, strong stellar lines 
in the vicinity of Ca I and Ca II K, as well as lower S/N, 
prevent any interstellar features from being discerned in 
this direction. Table 10 presents the equivalent widths 
of Ca I and Ca II K for the other two stars along with 
previous measurements and upper limits from the liter- 
ature. Figure 21 shows our fits to the observed Ca II K 
profiles. For a Leo, the fit results in a total column den- 
sity of A^(Ca II) = (0.69 ± 0.08) x 10^° cm'^. A similar 
fit to the Ca I profile yields a column density of N{Ca, I) 
= (0.42 ±0.11) X 10^ cm-2. For a Vir, our Ca II K profile 
synthesis fit gives iV(Ca II) = (7.11 ±0.15) x 10^° cm'^. 
No Ca I absorption is detected in the direction of a Vir, 
but we derive an upper limit on the column density of 
7V(Ca I) < 0.21 X 10^ cm-^. 
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Table 1 

Stellar and Observational Data 









B 


E{B - V) 




I 


6 




Exp. Time'' 




Star 


Name 


Sp. Type 


(mag) 


(mag) 


Ref. 


(deg) 


(deg) 


(pc) 


(hr) 


S/Nb 


McDonald Sight Lines 


HD 23408 


20 Tau 


B7 HI 


3.81 


0.07 


1 


166.17 


-23.51 


110 + 13 


15.5/14.5 


610/1240 


HD 23480 


23 Tau 


B6 IV 


4.11 


0.10 


1 


166.57 


-23.75 


110 + 13 


33.0/33.0 


1180/1930 


HD 24398 


C Per 


Bl lb 


2.97 


0.34 


2 


162.29 


-16.69 


301 + 68 


6.6/6.4 


810/1500 


HD 147933 


p Oph A 


B2 V 


5.26 


0.45 


3 


353.69 


+17.69 


90 


43.9/42.4 


1100/1600 


HD 149757 


C Oph 


09.5 V 


2.60 


0.32 


3 


6.28 


+23.59 


140 + 14 


7.5/7.5 


650/1590 


HD 179406 


20 Aql 


B3 V 


5.45 


0.33 


4 


28.23 


-8.31 


160 


22.3/22.8 


900/1230 


VLT/UVES Sight Lines 


HD 73882 


NX Vcl 


08.5 Vn 


7.58 


0.70 


4 


260.18 


+0.64 


820 


1.7 


1320 


HD 152236 


Sco 


Bl la+pe 


5.18 


0.68 


4 


343.03 


+0.87 


600 


0.2 


1090 


HD 154368 


V1074 Sco 


09.5 lab 


6.57 


0.78 


4 


349.97 


+3.22 


1100 


2.7 


1680 


HD 161056 




B1.5 V 


6.68 


0.59 


2 


18.67 


+11.58 


280 


3.1 


1480 


HD 169454 


V430 Set 


Bl la 


7.39 


1.08 


3 


17.54 


-0.67 


750 


6.4 


1990 


HD 170740 




B2 V 


5.93 


0.48 


4 


21.06 


-0.53 


213 + 43 


2.2 


1590 


HD 210121 




B9 V 


7.83 


0.31 


3 


56.88 


-44.46 


210 + 45 


1.3 


1160 



References. — (1) White et al. 2001; (2) Wegner 2003; (3) Valencic et al. 2004; (4) Rachford et al. 2009. 

Distances with error bars are from > 4-ct Hipparcos results. Other distances are derived by means of spectroscopic 
parallax. 

^ Exposure time and signal-to-noise ratio per resolution element for spectral region near CN A3874/CH"'" A4232. 
For the VLT/UVES sight lines, the values refer to the CN A3874 region. 



Table 2 

Comparison of Total Equivalent Widths for CH+ 
A4232 











Star 


(mA) 


(mA) 


Ref. 


20 Tau 


24.93 + 0.07 


0.36 + 0.07 


1 




23.3 + 0.3 


0.67 + 0.15 


2 




23.5 


0.57 + 0.11 


3 


23 Tau 


15.19 ± 0.05 


0.23 + 0.05 


1 




14.94 + 0.02 


0.15 


4 




14.0 + 0.2 


0.36 + 0.08 


2 


p Oph A 


12.50 + 0.06 


0.21 + 0.06 


1 




12.25 + 0.09 


0.11 + 0.05 


5 


COph 


23.00 + 0.06 


0.36 + 0.06 


1 




23.42 + 0.03 


0.29 


6 




23.8 + 0.09 


0.45 + 0.06 


7 




23.25 + 0.07 


0.41 + 0.02 


8 




22.62 + 0.04 


0.34 + 0.01 


9 




22.5 + 0.2 


0.54 + 0.1 


10 




21.7 + 0.2 


0.59 + 0.10 


11 




23.5 


0.36 + 0.06 


3 



References. — (1) This Work; (2) Hawkins & 
Jura 1987; (3) Vandcn Bout & Snoll 1980; (4) Stahl 
ct al. 2008; (5) Fcdcrman et al. 2003; (6) Casassus 
et al. 2005; (7) Hawkins et al. 1993; (8) Crane et 
al. 1991; (9) Stahl et al. 1989; (10) Hawkins et al. 
1989; (11) Hawkins et al. 1985. 
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Table 3 

Comparison of Total Equivalent Widths for CN A3874 







12CN 




13 


CN 






WxRlO) 




WxP{l) 


WxR{0) 


WxR{l) 




Star 


(mA) 


(mA) 


(mA) 


(mA) 


(mA) 


Ref. 


McDonald Sight Lines 


C Per 


8 88 + 06 


2 72 + 06 


1 36 + 06 


16 + 06 




1 




8.70 ± 0.6 


2.64 ± 0.2 


1.00 ± 0.3 






2 




9.06 ± 0.17 


2.79 ± 0.17 


1.36 ± 0.17 






3 




8.4 ± 0.21 






0.12 ± 0.02 




4 




8.99 ± 0.02 


2.89 ± 0.02 


1.30 ± 0.02 






5 




8.06 ± 0.07 


2.51 ± 0.04 


1.28 ± 0.04 






6 


p Oph A 


6.12 ± 0.04 


1.95 ± 0.04 


0.95 ± 0.04 


0.15 ± 0.04 




1 




6.53 ± 0.4 


2.18 ± 0.4 


1.33 ± 0.4 






2 




6.1 ± 0.4 


1.9 ± 0.3 


1.1 ± 0.3 






7 




6.11 ± 0.48 


2.08 ± 0.67 


1.05 ± 0.59 






8 




5.9 ± 0.3 


2.5 ± 0.3 


1.2 ± 0.7 






9 


COph 


7.98 ± 0.08 


2.37 ± 0.08 


1.27 ± 0.08 


0.19 ± 0.08 




1 




8.31 ± 0.3 


2.56 ± 0.4 


1.29 ± 0.6 






2 




8.17 ± 0.04 


2.50 ± 0.04 


1.29 ± 0.04 


0.27 ± 0.10 




3 




8.16 ± 0.03 






0.10 ± 0.05 




10 




7.75 ± 0.04 


2.45 ± 0.02 


1.26 ± 0.02 


0.19 ± 0.02 




11, 12 




7.34 ± 0.12 


2.13 ± 0.06 


1.08 ± 0.07 






6 




7.5 ± 0.3 


2.6 ± 0.3 


1.2 ± 0.7 






9 


20 Aql 


10.92 ± 0.07 


3.29 ± 0.07 


1.64 ± 0.07 


0.20 ± 0.07 




1 




11.16 it 0.2 


3.99 ± 0.2 


1.86 ± 0.2 






2 




11.5 ± 0.2 


3.2 ± 0.2 


1.6 ± 0.2 






13 


VLT/UVES Sight Lines 


HD 73882 


34.01 ± 0.06 


19.72 ± 0.06 


12.07 ± 0.06 


1.44 ± 0.06 


0.39 ± 0.06 


1 




34.83 ± 1.93 


17.31 ± 0.55 


10.23 ± 0.42 






8 




31.8 ± 0.7 


20.2 ± 1.3 


13.7 ± 1.2 


0.7 ± 0.5 




14 


HD 152236 (+6.0)'" 


7.09 ± 0.06 


2.22 ± 0.06 


1.20 ± 0.06 


0.13 ± 0.06 




1 


6.09 ± 0.32 


2.09 ± 0.60 


0.74 ± 0.59 






8 


HD 154368 (+5.2)=- 


25.30 ± 0.04 


15.59 ± 0.04 


9.66 ± 0.04 


1.13 ± 0.04 


0.33 ± 0.04 


1 


26.58 ± 0.6 


15.19 ± 0.7 


10.01 ± 0.2 






2 




25.34 ± 0.23 


14.85 ± 0.22 


9.21 ± 0.23 






3 




22.5 ± 2.8 


16.8 ± 2.3 


8.1 ± 0.9 






14 




24.62 ± 0.13 


15.12 ± 0.11 


9.74 ± 0.09 


0.82 ± 0.10 


0.25 ± 0.08 


15 


HD 161056 (-3.1) 


12.23 ± 0.05 


4.49 ± 0.05 


2.33 ± 0.05 


0.41 ± 0.05 




1 


11.47 ± 0.6 


3.38 ± 0.6 








2 




8.7 ± 1.0 


1.1 ± 1.0 


2.0 ± 1.0 






14 


HD 161056 (+2.8) 


14.15 ± 0.06 


4.26 ± 0.06 


2.22 ± 0.06 


0.43 ± 0.06 




1 


15.51 ± 0.6 


4.54 ± 0.6 








2 




15.6 ± 1.0 


4.7 ± 1.0 


3.0 ± 1.0 






14 


HD 169454 


22.22 ± 0.04 


16.50 ± 0.04 


11.51 ± 0.04 


1.92 ± 0.04 


0.47 ± 0.04 


1 




21.46 ± 0.2 


15.75 ± 0.3 


11.29 ± 0.2 






2 




25.10 ± 0.70 


15.50 ± 0.88 


10.85 ± 0.83 






8 




22.2 ± 1.0 


15.7 ± 0.4 


11.7 ± 0.5 


3.3 ± 1.3 


1.0 ± 0.5 


14 


HD 170740 


13.53 ± 0.05 


5.86 ± 0.05 


3.27 ± 0.05 


0.21 ± 0.05 




1 




13.60 ± 0.43 


6.68 ± 0.61 


3.82 ± 0.60 






8 




14.0 ± 0.4 


5.2 ± 0.8 


2.6 ± 0.3 






14 


HD 210121 


26.22 ± 0.07 


12.20 ± 0.07 


6.73 ± 0.07 


0.94 ± 0.07 


0.33 ± 0.07 


1 




27.9 ± 1.8 


10.5 ± 1.6 


4.5 ± 1.0 






16 




26.8 ± 3.0 


13.1 ± 2.6 


8.1 ± 2.4 






17 




23.4 ± 2.0 


10.0 ± 0.9 


4.9 ± 0.9 






14 



References. — (1) This Work; (2) Slyk et al. 2008; (3) Roth & Meyer 1995; (4) Kaiser et al. 1991; (5) Kaiser 
& Wright 1990; (6) Meyer & Jura 1985; (7) Pan et al. 2004; (8) Palazzi et al. 1992; (9) Federman et al. 1984; (10) 
Hawkins et al. 1993; (11) Crane et al. 1989; (12) Crane & Hegyi 1988; (13) Knauth et al. 2001; (14) Gredel et al. 
1991; (15) Palazzi et al. 1990; (16) Sheffer (private communication); (17) Welty & Fowler 1992. 

* Equivalent widths listed for HD 152236 and HD 154368 refer to the dominant velocity component. For the —4.8 
km s-^ component toward HD 152236, we find WxR{0) = 1.64±0.06, WxRW = 0.39±0.06, and W^P(l) = 0.14±0.06, 
while for the -13.9 km s~^ component toward HD 154368, we find WxR{0) = 0.28 ±0.04 and WxR{l) = 0.09 ±0.04. 
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Table 4 

Adopted Molecular Line Parameters 



Band 


Line 




/ 


Ref. 


CN B-X (0, 0) 


m 


3874.608 


0.0342 


1 




3873.998 


0.0228 


1 




p(i) 


3875.763 


0.0114 


1 




R(2) 


3873.369 


0.0205 


1 




P(2) 


3876.310 


0.0137 


1 


CN B-X (1, 0) 


R{0) 


3579.963 


0.0030 


1 


CH+ A-X (0, 0) 


m 


4232.548 


0.00545 


2 


CH+ A-X (1, 0) 


m 


3957.692 


0.00331 


2 



References. — (1) Roth & Meyer 1995; (2) Gredel et al. 1993. 



Table 5 

CN Equivalent Width Ratios and Derived Column Densities and 6- Values 



Star 


(0, 0) R(l)/P(l) 


jV(Af=l)^ 


b'' (0, 0)/(l, 0) iJ(0) 


N{N=0)'' 


6=- 


McDonald Sight Lines 


C Per 


1.98 ± 0.11 




oo 








p Oph A 


2.16 ± 0.10 




oo 


11.33 ± 3.36 




oo 


COph 


1.74 ± 0.17 




oo 


11.59 ± 1.54 




oo 


20 Aql 


2.01 ± 0.10 




oo 








VLT/UVES Sight Lines 


HD 73882 


1.64 ± 0.01 


9.84 ± 0.15 


0.94 ± 0.03 








HD 152236 (+6.0) 


1.91 ± 0.10 




oo 


13.28 ± 3.26 




oo 


HD 154368 (+5.2) 


1.62 ± 0.01 


7.91 ± 0.11 


0.72 ± 0.02 


4.65 ± 0.09 


18.60 ± 0.47 


0.68 ± 0.01 


HD 161056 (-3.1) 


1.86 ± 0.04 


1.67 ± 0.10 


0.52 ± 0.30 


15.64 ± 3.01 




oo 


HD 161056 (+2.8) 


2.00 ± 0.06 




oo 


16.57 ± 3.08 




oo 


HD 169454 


1.41 ± 0.01 


11.93 ± 0.14 


0.52 ± 0.01 


2.83 ± 0.05 


32.06 ± 0.88 


0.49 ± 0.01 


HD 170740 


1.77 ± 0.03 


2.42 ± 0.11 


0.42 ± 0.08 


7.18 ± 0.54 


5.86 ± 0.51 


0.47 ± 0.04 


HD 210121 


1.73 ± 0.02 


5.44 ± 0.16 


0.77 ± 0.08 


5.94 ± 0.16 


14.38 ± 0.46 


0.80 ± 0.02 



^ Column density (in 10^^ cm ^) and 6-value (in km s ^) derived from the measured equivalent width 
ratio for a pair of CN transitions arising from the same rotational state (either N=l or N=0). 
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Table 6 

CH+ and CN Component Structure 







CH+ 






CN 




Star 


(km s-1) 


TV* b 
(1012 cm-2) (km s-i) 


VLSR 
(km s-l) 


(1012 cm-2) 


(km s-l) 


McDonald Sight Lines 


20 Tau 


+7.4 


17.93 ± 0.09 


1.8 










+8.7 


17.20 ± 0.09 


1.8 








23 Tau 


+7.0 


1.06 ± 0.08 


2.7 










+10.0 


18.42 ± 0.06 


2.1 








C Per 


+6.7 


3.16 ± 0.08 


2.2 


+6.6 


2.31 ± 0.02 


0.93 


p Oph A 


+1.5 


6.95 ± 0.07 


1.8 


+1.7 


1.77 ± 0.01 


0.42 




+3.2 


8.89 ± 0.06 


1.5 








C Oph 


-0.9 


31.86 ± 0.08 


2.0 


-1.4 


1.61 ± 0.02 


0.60 










—0.2 


0.45 ± 0.02 


0.40 


20 Aql 


-6.6 


0.50 ± 0.09 


2.0 










+2.6 


3.15 ± 0.10 


2.4 


+1.7 


2.76 ± 0.02 


1.34 






VLT/UVES Sight 








LIU 1 OOO^ 


+5.6 


19.57 ± 0.09 


2.8 


-1-5 7 

-TO. / 


.uu HZ u.u^ 


0.88 




+10.2 


2.39 ± 0.07 


1.9 










Q 

— -Lu.y 


0.85 ± 0.15 


3.1 










-10.2 


3.38 ± 0.12 


2.0 










-2.0 


11.07 ± 0.15 


3.1 


-4.8 


0.38 lb 0.01 


0.54 




-4-7 1 


6.40 ± 0.14 


3.0 


+6.0 


9 91 _|_ n 09 


0.40 


HD 154368 


-13.0 


0.29 ± 0.06 


1.9 


-13.9 


0.06 ± 0.01 


0.20 




-6.8 


0.90 ± 0.06 


1.2 










+2.4 


14.44 ± 0.07 


2.1 










+6.6 


7.29 ± 0.07 


2.1 


+5.2 


18.19 ± 0.03 


0.68 ± 0.01 


HD 161056 


-7.9 


0.51 ± 0.15 


4.3 










-3.2 


0.77 ± 0.12 


3.5 


-3.1 


3.34 ± 0.01 


1.02 




+2.4 


24.31 ± 0.13 


3.8 


+2.8 


3.51 ± 0.02 


2.00 


HD 169454 


-3.0 


4.41 ± 0.08 


4.2 










+5.5 


14.96 ± 0.06 


2.9 


+5.7 


30.92 ± 0.05 


0.50 ± 0.01 




+15.4 


1.16 ± 0.05 


1.6 








HD 170740 


-1.1 


0.51 ± 0.06 


2.8 










+5.1 


16.46 ± 0.07 


3.1 


+6.6 


6.23 ± 0.02 


0.46 ± 0.04 


HD 210121 


-7.4 


2.82 ± 0.10 


1.5 


-6.7 


14.33 ± 0.04 


0.80 ± 0.02 




-2.7 


8.33 ± 0.13 


2.6 








^ Column densities refer to the "^^C-bearing isotopologue. 

^ We include 1-cr errors on b-values derived by means of the doublet ratio method (see 


§ 3). 
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l^ble 7 

Total Molecular Column Densities and ^^C/^^C Ratios 



Star 


Ar(i2cN) 
(1012 cm-2) 


12CN/13CN 


A?;(CN) 
(km s-l) 


Ar(i2CH+) 
(1012 cm-2) 


12CH+/13CH+ 


Ai;(CH+) 
(km s-i) 


Ar(i2CO) 
(10" cm-2) 


12CO/13CO 


Ref. 


McDonald Sight Lines 


20 Tau 


< 0.06 






35.13 ± 0.10 


82.5 ± 16.7 


-18.7 








23 Tau 


< 0.03 






19.48 ± 0.06 


72.7 ± 15.9 


-18.6 








C Per 


3.28 ± 0.03 


67.0 ± 28.0 


+13.2 


3.16 ± 0.08 






17.9 ± 0.5 


108 ± 5 


1 


p Oph A 
COph 


2.50 it 0.02 
2.89 ± 0.03=- 


54.0 ± 15.1 
48.8 ± 19.5 


+12.9 
+13.0 


15.84 ± 0.08 
31.86 ± 0.08 


64.7 ± 18.5 
75.2 ± 12.1 


-19.0 
-19.4 


19.2 ± 2.5 
25.4 ± 1.6 


125 ± 36 
167 ± 25 


1 
1 


20 A(il 


3.94 ± 0.04 


62.1 + 22.2 


+ 12.8 


3.64 + 0.13 






.30+^5 


■•50 ± 15 


2 


\'LT/i;VI£S Sight Lines 


HD 73882 


38.44 ± 0.06 


85.0 ± 3.3 


+12.6 


21.96 ± 0.11 






355 ± 208 


25 ± 22 


3 


HD 152236 
HD 154368 


3.52 ± 0.04^ 

26.98 ± 0.04 


64.9 ± 35.6 
70.7 ± 3.6 


+13.5'' 
+13.0 


21.70 ± 0.28 

22.91 ± 0.13 






26.7 ± 5.5 


37 ± 8 


1 


HD 161056 


10.10 ± 0.04 


36.3 ± 3.5 


+14.0 


25.59 ± 0.24 












HD 169454 


44.73 ± 0.06 


68.7 ± 1.3 


+12.8 


20.52 ± 0.12 












HD 170740 


8.77 ± 0.03 


133.6 ± 33.0 


+13.4 


16.97 ± 0.09 












HD 210121 


19.77 ± 0.05 


68.4 ± 4.9 


+13.3 


11.14 ± 0.16 






68 ± 14 


102 ± 31 


3 


Weighted mean 




67.5 ± 1.0 


+13.1 




74.4 ± 7.6 


-18.9 









References. — (1) Sheffer et al. 2007; (2) Hanson et al. 1992; (3) Sonnentrucker et al. 2007. 

" Total l^CN column densities for C, Oph and HD 152236 do not include contributions from tiie N = 2 level, but these contributions are 
likely to be small and within the uncertainties. If we assume T12 = Tfai, we get N{N = 2) = 0.02 x IO12 cm" ^ in both cases. 
The velocity shift for ON R{0) toward HD 152236 was held fixed during profile synthesis, and thus is not included in the mean. 



Table 8 

Comparison between Column Densities Derived from the (0, 0) and (1, 0) Bands 





Ar(CH+)a 


Ar(CN)=' 


Star 


(0, 0) R{0) 


(1, 0) R(0) 


(0, 0) R{0) 


(1, 0) R{0) 


McDonald Sight Lines 


20 I'au 


35.13 ± 0.10 








23 Tau 


19.48 ± 0.06 


20.16 ± 0.20 






C Per 


3.16 ± 0.08 




2.31 ± 0.02 




p Oph A 


15.84 ± 0.08 




1.77 ± 0.01 


1.53 ± 0.48 


COph 


31.86 ± 0.08 


33.24 ± 0.11 


2.06 ± 0.02 


2.13 ± 0.31 


20 Aql 


3.64 ± 0.13 




2.76 ± 0.02 




VLT/UVES Sight Lines 


HD 73882 


21.96 ± 0.11 


22.50 ± 0.26 


27.66 ± 0.04 




HD 152236 


21.70 ± 0.28 


21.42 ± 0.49 


2.59 ± 0.03 


2.36 ± 0.82 


HD 154368 


22.91 ± 0.13 


22.89 ± 0.24 


18.25 ± 0.04 


18.59 ± 0.38 


HD 161056 


25.59 ± 0.24 


27.05 ± 0.42 


6.85 ± 0.02 


5.44 ± 0.84 


HD 169454 


20.52 ± 0.12 


21.14 ± 0.20 


30.92 ± 0.05 


30.78 ± 0.52 


HD 170740 


16.97 ± 0.09 


17.62 ± 0.21 


6.23 ± 0.02 


5.64 ± 0.41 


HD 210121 


11.14 ± 0.16 


10.98 ± 0.40 


14.33 ± 0.04 


13.42 ± 0.38 



Total (line-of-sight) column densities (in I0I2 cm 2) derived through profile syn- 
thesis of the (0, 0) and (1, 0) R{0) lines in CH+ and CN. 



CN AND CH+ IN DIFFUSE MOLECULAR CLOUDS 



Table 9 

CN Rotational Column Densities and Excitation Temperatures 



Star 


N{N=0) 
(1012 cm-2) 


N{N=l) 
(1012 cm-2) 


N(N=2) 
(IQl^ cm-2) 


Toi 
(K) 


Tl2 

(K) 


roi(i3CN) 
(K) 


McDonald Sight Lines 


C Per 

p Oph A 
COph 
20 Aql 


2.31 ± 0.02 
1.77 ± 0.01 
2.06 ± 0.02 
2.76 ± 0.02 


0.94 ± 0.02 

0.68 ± 0.01 
0.82 ± 0.02 
1.12 ± 0.02 


0.04 ± 0.02 
0.04 ± 0.01 

0.06 ± 0.02 


2.723 ± 0.031 

2.657 ± 0.026 
2.702 ± 0.042 
2.728 ± 0.028 


2.931 ± 0.411 

3.136 ± 0.323 

3.155 ± 0.340 








VLT/UVES Sight L 


ncs 






HD 73882 
HD 152236 (-4.8) 
HD 152236 (+6.0) 
HD 154368 (+5.2) 
HD 161056 (-3.1) 
HD 161056 (+2.8) 
HD 169454 
HD 170740 
HD 210121 


27.66 ± 0.04 
0.38 ± 0.01 
2.21 ± 0.02 

18.19 ± 0.03 
3.34 ± 0.01 
3.51 ± 0.02 

30.92 ± 0.05 
6.23 ± 0.02 

14.33 ± 0.04 


10.48 ± 0.02 
0.12 ± 0.02 
0.81 ± 0.02 
8.40 ± 0.02 
1.59 ± 0.02 
1.46 ± 0.02 

13.30 ± 0.02 
2.44 ± 0.02 
5.31 ± 0.02 


0.31 ± 0.02 

0.29 ± 0.01 
0.10 ± 0.02 
0.10 ± 0.02 
0.51 ± 0.01 
0.11 ± 0.02 
0.13 ± 0.02 


2.631 ± 0.004 
2.454 ± 0.178 
2.588 ± 0.033 
2.911 ± 0.004 
2.960 ± 0.018 
2.761 ± 0.020 
2.804 ± 0.003 
2.672 ± 0.011 
2.604 ± 0.007 


2.693 ± 0.038 

2.815 ± 0.033 
3.333 ± 0.159 
3.383 ± 0.210 
2.881 ± 0.017 
2.989 ± 0.123 
2.589 ± 0.096 


2.706 ± 0.200 

2.781 ± 0.192 

2.599 ± 0.082 
3.096 ± 0.383 


Weighted mean 








2.754 ± 0.002 


2.847 ± 0.014 


2.652 ± 0.069 



Table 10 

Total Ca I and Ca II Equivalent Widths 



Star 


WxiCa. I) 
(mA) 


WxiCe. II) 
(mA) 


Ref. 


a Leo 


0.12 ± 0.03 


0.59 ± 0.07 


1 






^ 0.4 


2 


a Vir 


< 0.06 


5.95 ± 0.12 


1 




< 0.2 




3 






5.4 


4 






4.8 ± 0.1 


2 






5.2 


5 






5.7 


6 






5 


7 



References. — (1) This Work; (2) 
Vallorga ct al. 1993; (3) Wclty ct al. 
2003; (4) Welty et al. 1996; (5) Lalle- 
ment et al. 1986; (6) Hobbs 1978; (7) 
Hobbs 1975. 



